TEW #7 g-H-xWH 


ANALYSIS OF AIRCRAFT WING-MOUNTED ANTENNA PATTERNS 

R. J. Marhefka 



TECHNICAL REPORT 2902-25 
June 1976 


Grant NGL 36-008-138 


The material contained in this report is also used as a dissertation 
submitted to the Ohio State University Electrical Engineering 
Department as partial fulfillment of the requirements for the degree 
Doctor of Philosophy. 


National Aeronautics and Space Administration 
Langley Research Center 
Hampton, Virginia 23665 



ACKNOWLEDGMENT 


The author would like to express his sincere gratitude to 
his graduate adviser. Professor R. C. Rudduck, for all the 
guidance and encouragement that he has provided. The assistance 
and cooperation of all his associates at the ElectroScience 
Laboratory are greatly appreciated. Deserving special recognition 
is Dr. W. D. Burnside for all his stimulating discussions and 
suggestions. Thanks must, also, be extended to Dr. P. H. Pathak 
and Dr. C. L. Yu for their help and ideas. Professor C. H. 

Walter and Professor L. Peters have been most helpful in their 
critical review of the author's manuscript. 

The efforts of Messrs. G. Ryno and J. Huang of the Naval 
Weapons Center, China Lake, California; and Messrs. W. Croswell 
and M. Gilreath of the National Aeronautics and Space Admin- 
istration, Langley Research Center in providing measured results 
on scale model aircraft are greatly appreciated. 



ABSTRACT 


The analysis of high frequency radiation patterns of aircraft 
wing mounted antennas is the object of this research. This is an 
analytical study of basic antenna types using ray optical techniques. 

The aircraft is modelled in its most basic form so that this study 
is applicable to general -type aircraft. The fuselage is modelled as 
a perfectly conducting finite elliptic cylinder. The wings and 
horizontal and vertical stabilizers are modelled as perfectly con- 
ducting "n" sided flat plates that can be arbitrarily attached to 
the fuselage or to themselves. Presently, the antenna locations are 
assumed to be on the surfaces of the wings at locations removed 
from engines and stores such that these effects are negligible. 
Volumetric patterns are calculated for several aircraft. The validity 
of the solution is shown by comparing the results against scale model 
measurements. The application of this solution to practical airborne 
antenna problems has shown its versatility in designing antennas and 
predicting their radiation patterns in an accurate and efficient manner. 
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CHAPTER I 
INTRODUCTION 


In order for modern airborne antenna systems to function prop- 
erly, the antenna pattern must meet certain specifications. In 
fact, system performance often is very much dependent on the 
pattern performance of the antenna. This requires the aircraft 
antenna designer to not only concern himself with the development 
of an antenna element for a particular application, but to also be 
concerned with the performance of the antenna system in the presence 
of the aircraft structure. The location of the antenna on the air- 
craft, and even the element design, is critical if a reliable 
antenna system is to be achieved. An efficient and economical 
means of predicting the pattern performance of an antenna is, there- 
fore, desirable. 

The object of this research is to analyze radiation patterns 
for wing mounted antennas at high frequency. Wing locations have 
become increasingly desirable for many applications, as the avail- 
able real estate for antennas becomes more difficult to find because 
of the large number of antenna systems on a modern aircraft. This, 
then, is an analytical study of basic antenna types for frequencies 
above 100 megahertz in which the antenna is mounted on or near the 
wing. 


The aircraft is modelled in its most basic form so that this 
study is applicable to general-type aircraft. The fuselage is 
modelled as a perfectly conducting finite elliptic cylinder. The 
wings and horizontal and vertical stabilizers are modelled as 
perfectly conducting "n" sided flat plates that can be arbitrarily 
attached to the fuselage or to themselves. Presently, the antenna 
locations are assumed to be on the upper surfaces of the wings or 
at locations removed from engines and stores such that these effects 
are negligible. The theory, however, is applicable to the addition 
of engine and store effects by modelling these structures as flat 
plates or small finite elliptic cylinders. 

This type of analytic solution enables antenna designers to 
investigate an antenna system's performance in the presence of the 
aircraft structure using a computer simulation method, in a manner 
that can save a great deal of engineering time and money. This 
allows airborne antennas to be located on the aircraft at the design 
stage of the aircraft. This means that optimum locations and nec- 
essary structural changes can be anticipated. Also, a future re- 
location or addition of antennas on an aircraft within its useful 
lifetime can be easily accomplished. Once an optimum region is 
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determined using computer simulation models of the aircraft, the 
antenna can be experimentally tested to ascertain its actual per- 
formance. 

It is not uncommon for antennas to be allocated specific 
locations which are not related to electronic system requirements. 
Typical ly, these locations are chosen primarily for convenience 
with regard to aircraft structural specifications. Scale model 
measurements, then, are used to evaluate the performance of 
the antenna system in terms of its desired pattern. This approach 
can be very time consuming and costly. 

The computer simulation method can, also, be used to determine 
optimum antenna designs for a given application. Various antenna 
types can be easily examined on the analytical models and the 
results compared. These analyses are based on arbitrarily- 
oriented electric and magnetic delta function sources which 
can be used to solve for the pattern of an arbitrary antenna 
simply by integrating over the equivalent aperture currents. 

Three basic methods have been used to compute on^aircraft 
antenna patterns. Among the first methods used were the modal 
solutions for infinitely long circular [1,2] and elliptic [3] 
cylinders. These solutions modelled the fuselage by a cylinder 
whose cross-section approximated the fuselage cross-section at the 
antenna location. Arbitrary antennas were considered in these 
studies in which the antenna was mounted either on or over the 
fuselage. This method could not consider the effects of the other 
scattering structures of the aircraft such as the wings and 
stabilizers. The effects of these additional structures can be 
quite large, in general. 

The second method for analysis of aircraft radiation patterns 
is primarily a low frequency technique. This approach uses a 
moment method technique to solve for the unknowns to an integral 
equation. The surface currents and resulting scattered fields can 
be found by enforcing the boundary conditions on the aircraft struc- 
tures. One of the first moment solutions applied to aircraft prob- 
lems was the wire grid technique, developed by Richmond [4], which 
employed a point matching technique [5]. This solution required the 
determination of approximately 100 unknown currents per square wave- 
length in order that the wire grid adequately model the perfectly 
conducting surface. Richmond [6] has developed a more sophisticated 
approach in which the reaction technique is used to solve for the 
unknown currents. This solution still requires approximately 100 
unknown currents per square wavelength. Another approach is to 
divide the surface of the conducting body into patches each having 
two orthogonal unknown complex currents [7]. This surfaced-current 
approach reduces the unknown currents to about 20 per square wave- 
length. This allows a larger surface to be considered. However, 
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all these solutions are restricted to low frequencies due to 
the limitation on the size of matrices which modern computers 
can solve without excessive loss of accuracy. 

The approach that has been taken in this study is the 
Geometrical Theory of Diffraction (GTD). This is basically 
a high frequency method that is applied in terms of two character- 
istic problems; these being wedge diffraction and curved surface 
diffraction. The wedge diffraction solution has been applied 
to determine the radiation patterns of such basic antennas as 
parallel plate antennas [8,9,10], parallel plate arrays [11,12], 
horn antennas [13,14], parabolic reflectors [15,16], and rectangular 
waveguide antennas [17]. Both of these types of diffraction 
solutions have been applied in computing the patterns of antennas 
mounted on cylinders [18,19,20] and rockets [21]. GTD has also 
been successfully applied to a number of problems involving 
antennas mounted on the fuselage of aircraft [22<-30], Antennas 
mounted on various isolated wing structures have been analyzed 
[31,32], however, these did not include the effects of the fuselage 
and other structures. 

Using this approach, one applies a ray optics technique 
to determine components of the field incident on various scatterers. 
Components of the diffracted field are found using the GTD solutions 
in terms of rays which are summed with the geometrical optics 
terms in the far field. The rays from a given scatterer tend 
to interact with the other structures causing various higher 
order terms. In this way one can trace out the various possible 
combinations of rays that interact between scatterers and determine 
and include only the dominant terms. Thus, one need only be 
concerned with the important scattering components and neglect 
all other higher-order terms. This makes the GTD approach ideal 
for a general high frequency study of aircraft antennas in that 
only the most basic structural features of the aircraft need 
to be modelled. - 

Since GTD is essentially a high frequency solution, the 
lower frequency limit of this solution is dictated by the spacings 
between the various scattering centers. In practice, this means 
that the scattering centers should be greater than approximately 
a half-wavelength apart. In some cases even this requirement 
can be relaxed. For a general type aircraft, therefore, the 
low frequency limit will be around 100 MHz. The upper frequency 
limit is dependent on how well the theoretical model simulates 
the important details of the actual structure. 

The basic approach applied in this study of antennas mounted 
on the wings of aircraft is to separately analyze the basic 
scattering shapes making up the computer simulated model of 
an aircraft. The single structures are analyzed using the ray 
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optics approach such that only the dominant effects are included. 

The scattered fields of the various component parts of the structure 
are found using the Geometrical Theory of Diffraction. The 
basic solutions needed for this study are briefly discussed 
in Chapter II. This includes geometrical optics reflected fields 
from curved surfaces, the diffraction by a wedge and curved 
wedge, and the slope diffraction by a wedge using arbitrarily 
oriented dipole moment sources. A heuristic corner diffraction 
coefficient is also introduced. Curved surface diffraction 
is presented along with a newly developed transition function 
for the shadow boundaries of curved surfaces. 

Chapter III describes the scattering from a number of finite 
flat plates. They can have an arbitrary orientation in space 
and can be attached to form bent plates. These plates are used 
to model the wings, horizontal and vertical stabilizers. 

Chapter IV describes the scattering from a finite elliptic 
cylinder. This is used to model the fuselage and can be used 
to model engines and stores. 

The complete aircraft model is described in Chapter V. 

The dominant scattering terms due to the interaction of the 
plates and cylinder are discussed. The analytic results are 
compared against scale model measurements to determine the validity 
of the approximations made in the analysis. Finally, a summary 
of the present study and a discussion of future topics are made 
in Chapter VI. 
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CHAPTER II 

THEORETICAL BACKGROUND 


A. Introduction 

The analytical modelling of an aircraft to predict the 
radiation patterns of antennas mounted on the wings has been 
accomplished by the use of the Geometrical Theory of Diffraction 
(GTD). This is a high frequency technique that allows a complicated 
structure to be approximated by basic shapes representing canonical 
problems in the GTD. These shapes include flat and curved wedges 
and convex curved surfaces. The GTD is a ray optical technique 
and it therefore allows one to gain some physical insight into 
the various scattering and diffraction mechanisms involved. 
Consequently, one is able to quickly seek out the dominant or 
significant scattering and diffraction mechanism for a given 
geometrical configuration. This, in turn, leads to an accurate 
engineering solution to practical antenna problems. The basic 
GTD solutions needed are discussed in this chapter. These solutions 
are applied to specific structural scattering problems in the 
following chapters. 

Two of the basic diffraction problems are illustrated in 
Figures la and b. In the format of GTD, the total electric 
field may be represented as 

E = EV + E r u r + E* 1 . 

The field E 1 is the electric field of the source in the absence 
of the surfaces. The field E** ts the electric field reflected 
from the surfaces. The field E^ is the diffracted field from 
the discontinuities of the structure, such as edges, corners, 
curved surfaces, etc. The step functions u 1 and ur are shown 
to emphasize the discontinuity in the incident and reflected 
fields at the incident and reflection shadow boundaries. They 
are one in their lit regions and zero in their shadow regions. 

The extent of the regions are determined by geometrical optics. 

The surfaces used in this report are perfectly conducting and 

the surrounding medium is free space. An exp (jut) time dependence 

is assumed throughout and suppressed. 
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OBSERVATION 


LIT REGION 



BOUNDARY 


(a) BASIC WEDGE PROBLEM 



Figure 1--Basic diffraction problems. 
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B. Geometrical -Optics Fields 


The incident electric field E 1 can be produced by an arbitrary 
electric or magnetic source causing plane, cylindrical, conical, or 
spherical wave incidence on the scattering structures. In this 
report the spherical wave case is treated. The far field pattern of 
the source can be written as 


E 1 ^,^) ~ C Q F (9 »4> ) + <l>G(e ,<}>)] 


e~ jks 

s 


( 1 ) 


The reflected electric field from a curved surface S as shown in 
Figure 2 is given in terms of geometrical optics by 


E^s) = E 1 (Q r ) 


r r 
p l p 2 


(p^+s)(p£+s) 


-jks 


( 2 ) 


where E^C Qr ) is the incident electric field at the point of reflection 
Q r and R is the dyadic reflection coefficient such that 


e,. - 


e x e x 


(3) 


The vector e^ is th§ unit vector perpendicular to the plane of in- 
cidence and §,i and e£ are the unit vectors parallel to the plane 
of incidence. The point of reflection Qr is found from the laws of 
reflection that state that the angle of incidence is equal to the 
angle of reflection that is 

-I • n = s • n (4a) 

and that the incident, reflected, and surface normal vectors must lie 
in the same plane or in other words 

A A A A 

I x n = s x n. (4b) 

r r 

The quantities p-j and p 2 are the principal radii of curvature of the 
reflected wavefront at the reflection point Qr. Kouyoumjian [34] 
discusses how to find these values for an arbitrary wavefront by 
diagonalizing the curvature matrix for the reflected wavefront 
given by Deschamps [35], Let a wavefront be incident on a curved 
surface S at Qr as shown in Figure 2. ei and e2 are unit vectors 
in the principal directions of.S at Qr with principal radii of 
curvature R-j and R 2 . X^j and are the principal directions of 
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in which the plus sign is associated with pT and the minus sign with 
p£ and where 



(6a) 


| 0 |= 9^022 - e 12 9 21 * 


(6b) 


The principal directions of the principal radii of curvature of the 
reflected wavefront are given by 

"<£,-^1 £-<£, 


<22 r A 1 w 12 2 
P 1 


X r = x 


x£ = -s x 


(7a) 


(7b) 


where 


•A A , 


= X n ' 9 - 2(n • X-j ^)n 


' 1,2 1 ,2 


(8) 


and 


<11 i 


1 , 2 cose 1 

- 2 
6 


22 


21 


2 J 


(9a) 


>r _ -2 cose 1 


<12 


0 22 9 1 2 + 9 11 0 21 


(9b) 


n r _ 1 , 2 cos 0 

^22 " i + i „ , 2 


p 2 


er L"! 


’12 


) Z 1 

hi 


(9c) 
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If the incident wavefront is spherical, the equations for the 
principal radii of curvature of the reflected wavefront are 
gi ven by 





O 

O 1 


1 

. JL + . 

1 

sin 0 £ sin^O-j 


r 

P 1 ,2 

" s' + 

cose 1 

L R ! 

R 2 J 



1 

' . 2 
sin 0 

. 2 

2 sin e.| 

2 

I f 

1 

- 

\ 

cos 2 © 1 ' 

R i 

+ R 2 J 

1 R 1 R 2 


where s' is the radius of curvature of the incident wavefront at Qp, 

0] is the angle between the direction of the incident ray I and unit 
vector e-| and 02 is the angle between the direction of the incident 
ray 1 and unit vector e^. 

C. Diffraction by a Wedge 

The wings, vertical and horizontal stabilizers on an aircraft can 
be approximated as flat plates. The Geometrical Theory of Diffraction 
provides a means of including the three-dimensional effects of the flat 
plates with straight edges by the use of wedge diffraction. This 
section briefly discusses this basic diffraction problem. 

An asymptotic solution for the diffraction from a conducting 
wedge was first solved by Sommerfeld [36]. Originally, GTD 
[37] as applied to diffraction by a wedge was based on plane 
wave diffraction coefficients; however, as shown in Reference [38] 
the use of cylindrical wave diffraction coefficients has been found 
necessary in the treatment of antennas. Consequently, different 
formulations of wedge diffraction were substituted for the plane 
wave diffraction coefficient which is the basis for wedge diffraction 
theory. Pauli [39] introduced the Vb function as a practical formu- 
lation to the solution for a finite-angle conducting wedge. Hutchins 
and Kouyoumjian [40,41], however, have presented a formula for 
the diffracted field, which significantly improves the accuracy 
over that obtained from Pauli's form. 

This improved diffraction solution [40,41]' is better in the 
transition regions (near the incident and reflected shadow boundaries). 
It can be written in the form 



where 


, -j (kL+ir/4) / \ 

! ±„< l - 6 - n) ' — 7r= J 7 cos (“if) 

jny 2 tt \ / 


x e 


jkLa 


• i- 

e" JT dx + [higher order terms] 


kLa 


and where the higher order terms are negligible for large kL and 
with n defined from the wedge angle WA = (2-n)ir, also 

a = a 1 = 1 + cos^-ZnirN 1 ) and 


N + is a positive or negative integer or zero, whichever most nearly 
satisfies the equations 

2mrN - 6 = - it for I 

-7T 

2mrN + - B = + it for I . . 

• 7T 

The variables L and 8 are defined below. 

The three-dimensional wedge diffraction problem is pictured in 
Figure 3. A source whose radiated E field is given by P(s) is 
located at point s' (p 1 ,<t»' ,z' ). It can be an arbitrary electric or 
magnetic source causing plane, cylindrical, conical, or spherical 
wave incidence on the wedge top. The diffracted vector field at 
point s(p,<j>,z) can be written in terms of a dyadic diffraction coef- 
ficient. Kouyoumjian and Pathak [42] have given a more rigorous 
basis for the GTD formulation and have shown that the diffracted 
fields may be written compactly if they are in terms of a ray-fixed 
coordinate system. The ray-fixed coordinate system is centered at 
the point of diffraction Qe, (or points of diffraction in the case 
of plane wave incidence). Qe is a unique point or points for a given 
source and observation point. The incident ray diffracts as a cone of 
rays such that e Q = (see Figure 3). 

The relationships between ^the orthogonal unit vectors associated 
with these coordinates (s 1 ,8^,^' ;s ,B 0 ,^) are given by 

I = - s' 

I = Bq x <Jj' 

A A A 

s = e 0 x <(., 


n 




where I is the incident direction unit vector, and s is the dif- 
fraction direction unit vector. The diffracted field is now given 
by 

E^s) * E i (Q E ) • f E (s,i)A(s)e" jks (11) 


where 


and 


D r = -$'8 D - 4>'<j>D, 
E ^o o s * * h 


_p-jV4 

D s h ( W'.?o) = -7i: . 

s,n 0 sine L 


cot P^- F[kLa+ (8 - )] 


+ cot^-j F[kLa-(e")] 


+^cotft^-^ F[kLa + (6 + )] + cot[^i^— ^ F [ kLa ” ( 8 + )] ^ 


2n 


J 


( 12 ) 


v. 


and 


F(X) = 2j|J)f|e jX | e~ jT ^ dx 

h[xl 


(13) 


is called the transition function. This can be shown to be related 
to the Vg function for a straight wedge as 


, rr _ jkL 

D Sjh U,V,B 0 ) = CV B (L.S') ♦ V B (L,B + )]^y— 
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In matrix notation this can be written as 


'Ej(s) 

_eJ(s) 




A(s)e 


-jks 


The Ds coefficient (minus sign between Vg terms) applies for the 
E-field component parallel to the edge with the boundary condition 
(acoustically soft) 


I wedge ^ ~ 


The Dh coefficient (plus sign between Vb terms) applies for the E- 
field vector perpendicular to the edge with the boundary conditions 
(acoustically hard) 



The angular relations are expressed by 

b = e + = <i>+ <f>' . 


The terms are associated with the incident field and the <j>+<j>' 
terms are associated with the reflected field. The quantity A(s) 
is a ray divergent factor given in general by [34] 

• Jiukr 

where for a wedge p = pi which is the radius of curvature of the 
incident wavefront in tne plane containing the incident ray and the 
edge. For ^ome specific problems this is given by 

plane, cylindrical (s=p), and conical 
wave incidence 


spherical wave incidence. 


A(s) - 


1 
1/7 


s’ 


s(s *+s) 
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The quantity L is a distance parameter for a straight wedge given 
in general for a straight wedge by [34] 


• • • Q 

S ^ p e +S ) p i p 2 sin B o 

pg(p]+s)(p2+s) 


(14) 


where p] , p \ are the principal radii of curvature of the incident 

wavefront. For some specific problems this is given by 
r 

2 

s sin plane wave incidence 





p 1 p 

p+p' 

s's sin 2 B 
. o 

s+s ' 


cylindrical wave incidence 


conical and spherical wave 
incidence. 


At grazing incidence, that is when <j>' = 0, the expression for 
Dfo and D s must be multiplied by a factor of 1/2, this comes about 
because the incident and reflected field merge together and only 
one-half the total field on the surface is the incident field, the 
other half is associated with the reflected field. 

D. Slope Diffraction by a Wedge 

The diffraction coefficient discussed in Section C is only 
valid for incident fields on an edge that are slowly varying in 
amplitude. If this is not true, as when a dipole is aligned such 
that its axis points in the direction of the edge, an additional 
term is needed that relates the diffracted field to the spatial 
derivative of the incident field. This term is referred to as 
the slope diffracted field. 

Rudduck and Wu [43] have previously obtained a slope diffracted 
solution (Vs) for the two dimensional general wedge for the hard 
boundary condition by differentiating the eigenfunction solution. 
Karp and Keller [44] first introduced the GTD slope diffracted 
field; however, this result is not valid at or near the shadow 
boundary. Recently Mentzer, Peters, and Rudduck [45] obtained 
a two-dimensional slope diffraction function for the special case 
of the half plane by differentiating the exact solution given 
by Sommerfeld. They applied this to horn antennas. Kouyoumjian 
and Hwang [46,47] have extended the results of Reference [44] 
such that Equation (15) is now valid even within the transition 


15 

I 


s 



region adjacent to shadow boundaries. The result presented here 
is for the diffraction coefficient for the straight wedge case 
given in Equation (12). 

This dyadic slope diffracted field for a dipole whose axis is 
directed radially toward the edge, that is, such that the amplitude 
of the incident field is zero at the edge but with a finite incident 
slope is given by [46,47] 

rd . J- • ^Afsle-^ (15) 

an jk 3$' 


where 


aD F „ , 3D A * 3D. 

a^ = - Wo W - 


The slope diffraction coefficient is explicitly given by 

csc^W kLaV,) 


-i - 

aD . J 4 

s,h = _ e 


a<j>' 


4n 2 N f27ik sin 2 g 


V 


- csc2 ^-j F s (kLa"(g")) 


|csc 2 (^j F s (kLa + (B + )) - csc 2 ^) F s (kLa " (3+)) 
where B + and a + are defined as in Section C and 


F S (X) 


= 2jX 


+ 4X 3/2 e :iX 




i/ri 


or 


F $ (X) = 2jX[l - F ( X ) ] 
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with F(X) defined as in Equation (13). 

The partial derivative of the incident field (a/an) is defined as 
the derivative taken in the direction normal to the edge-fixed plane 
of incidence. A very important source type to be considered for 
general application is the arbitrarily oriented electric and magnetic 
dipole moment which can be expressed by 

u'PgSds-s* |). 
m 

Any general source can be approximated by superimposing a weighted 
set of these dipole moment elements to obtain a sufficiently accurate 
result for engineering purposes [48]. 

If the dipole moment has an arbitrary orientation in space with 
respect to the wedge, the ordinary diffraction effects and the slope 
diffraction effects should both be considered for completeness. In 
matrix notation, this may be written, assuming a spherical wavefront 
incidence, as 


En(s) 

A . 

1 aD S 3 

_n - S3 q 

s jk a* 1 an u 


e1«e> 

Ei(s)_ 

* U 

i 3D. a 

n _n JL \L iL_ 

L h jk a*' an 


e 1 ( q e ) 


x 


s 1 



(16) 


The dipole moment producing the incident field can be decomposed 
into three components in the ray-fixed coordinate system. Thus 


p e = p e i + p e 




m 


ml 


mg' 

o 


m 


<f>' 




where 




Pe 

m 
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I 



= B, 


Pe 

m 




e . 
m 


The electric field of an electric dipole moment (TM case) can be 
written in the far field as 


Eg = njk p e sin e e 
P 


-jks 

4its 


where n is the 
electric field 
as 



impedance of surrounding media. The far field 
of a magnetic dipole moment (TE case) can be written 


m 


sin e p 


£jks 

4irS 


where 0n, 4 > p are measured with respect to the dipole coordinate sys- 
tem. Tne non-zero incident fields needed to compute the diffracted 
field for the electric dipole moment case are 


Em = B' • E 1 = njk p Q0 , 


-jks 


eg' 


4irs 1 


E 1 = i' ’ E 1 = njk p , 


-jks 1 


e4> ' 4 its 1 


an 


aE 


= r 


an 


= njk p 


-jks 1 


e I 4 tts' 2 


The non-zero incident fields needed to compute the diffracted fields 
for the magnetic dipole moment case are 
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g-jks 1 
4tt S ' 


El 



p m0' 


e" J ’ ks ' 
4tts 1 


^E|", a gpi _ p“jkS 

= 6 ' ' ”aiT -jk p ml ^X 


The above fields may be used in Equation (16) to give the dif- 
fracted fields from the edge. 

E. Diffraction by a Curved Wedge 

The fuselage of a general aircraft is to be modelled in this 
report as a finite elliptic cylinder. The ends of the cylinder 
are formed by a flat surface and the elliptic cylinder surface 
with a resulting curved edge at their junction. The diffracted 
field due to this curved wedge is analyzed using the GTD techniques 
developed by Kouyoumjian and Pathak [34]. 

The curved wedge problem is illustrated in Figure 4. The dif- 
fracted field from the curved wedge may be written in the form 


^(s) 


E 1 


(Q e ) 


D e (s,I) 


J 


s(p+s) 


,-jks 


(17) 


The parameter p is the distance between the caustic at the edge and 
the second caustic of the diffracted ray. This is given by 


1 = 1 _ 

p pi 


n e • (I - i) 
a e s, ' nife o 


(18) 


where Pg, n e , a e are defined below. The diffraction coefficients 
for the curved wedge are extended from those in Equation (12) to 
allow the diffracted field to be continuous at the incident and 
reflected shadow boundaries. This is accomplished by finding the 
appropriate distance parameters L in each of the transition functions 
that make the field continuous. The diffraction coefficient for the 
curved wedge is given by 
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Figure 4--Geometry for three-dimensional curved 
wedge diffraction problem. 
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(+,+', 8 0 ) 


e~ j7T/4 

2n \/2irk sing Q 


‘2 sin J F[kL 1 a(B"}] 

- cos (w) - cos (\) 



(19) 


in which a(g) = 2 cos 2 6/2 and a + (g) = 2 cos 2 . The 

distance parameter associated with the incident field is given by 

- 2 's 

. . . (20a) 

p e( p l +S ) ( p 2 +S ) 


L i = s ( p e +s )Pl p 2 siA o 


The distance parameter associated with the reflected field from the 
surface with superscript 0 (reflection boundary at tt-V) and the re- 
flected field from the surface with superscript n (reflection bound- 
ary ( 2n-l )tt-«j> * ) is given by 



(20b) 


The parameters p] and p 2 are the principal radii of curvature of 
the incident wavefront at the diffraction point Qe- The parameter 
p^ is the radius of curvature of the incident^wavefront at Qe taken 
in the plane containing the incident ray and e the unit vector 
tangent to the edge at Qe- In the case of plane, cylindrical, and 
conical waves pi is infinite and for the case of spherical waves 
pi=s'. The parameters p^ and pj? are the principal radii of 
curvature of the reflected wavefront at Qe> which may be found from 
Equation (10) for spherical wavefront incidence. The parameter p£ 
is the radius of curvature of the reflected wavefront at Qe taken 
in the plane containing the reflected ray and e the unit vector 
tangent to the edge at Qe* This can be found from the expression 


21 



1 1 


where n is the unit normal to the surface at Q^, n e is the associated 
unit normal to the edge directed away from the center of curvature, 
and a g >0 is the radius of curvature of the edge at Q^. 

A slope wave term can be found for the curved wedge problem; 
however, for the geometries involved with wing mounted antennas this 
term does not appear to be a significant contributor to the total 
scattered field so it will not be used in this study. 

F. Diffraction by a Corner 

The wings of aircraft are modelled in this study as flat plates. 
These plates can be thought of as being composed of a series of finite 
straight edges that join together at the corners of the plate. The 
diffraction of energy from these edges, as discussed in Section C, 
is based on diffraction from infinite straight edges. To com- 
pensate for the finiteness of the edges, a diffraction coefficient 
associated with the corners of the plates is needed. 

An eigenfunction solution for the scattering from the corner of 
wedges has been derived by Satterwhite [49]. This solution is ex- 
tremely cumbersome for numerical calculations. A GTD diffraction 
coefficient is needed that is numerically efficient in order for 
the corner effect to be of any practical use in complicated 
modelling problems. Burnside and Pathak [50] have recently proposed 
a solution based on the asymptotic evaluation of the radiation in- 
tegral which employs the equivalent edge currents that would exist 
in the absence of the corners. The corner diffraction term is then 
found by appropriately (but presently heuristically) modifying the 
asymptotic result for the radiation integral which is characterized 
by a saddle point near an end point. This diffraction coefficient 
is still in the beginning stages of its development. However, it 
has been shown to be useful in eliminating discontinuities in the 
total field associated with the corner effect. For this reason, 
it is briefly discussed here and has been included in the results 
as a good engineering approximation to the problem. 

The corner diffraction problem is illustrated in Figure 5. The 
fields associated with one edge of the corner diffraction fields in 
the far field with spherical wave incidence are given by 


2(n-n e )(I*n) 


a e sin2|S o 


(21) 
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(22a) 
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dnd 
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C s,h (Q E ) 


. 7T 

-J 4 


2j2irk sing ( 


r 


r ~\ 


. F[kLa(3‘) 

p 

a(B') 1 


COS 

^ L ' 

r 

_ks c a{6o-S c ) 1 



F[kLa(3 )] 

+ + 

cos ^2 


a(e + ) 


. ks c a( V B c } -i 


(22c) 


The function F(x) is defined in Section C by Equation (13), and 
as before a(3) = 2 cos^ 3/2 where 3 ±= T 4> ' » and L = s' sin^ 3 0 for 
spherical wave incidence. The function C s ^(Qe) 1S a modified 
version of the diffraction coefficient in liquation (12) for the half- 
plane case when n = 2. The modification factor. 


a( 

,e) t 

ks a 
c 



F(kLa), 


is a heuristic function that insures that the diffraction coefficient 
will not change sign erroneously when it passes through the shadow 
boundaries of the edge. The F(kLa) in the I and M terms is also 
heuristic in that it is this term along with the function 



that gives the sign change. There is also a corner diffraction term 
associated with the other edge forming the corner and is found in 
a similar manner. 

Even though further study is needed to improve the diffraction 
mechanism, it is felt that the benefits obtained from its inclusion 
warrants its use here. Experimental results in Chapter III, Section 
B, confirm its utility in engineering applications. 

G. Curved Surface Diffraction 

The aircraft fuselage is a long cylindrical type structure with 
a generally convex curved surface. The GTD mechanism that describes 
the scattering in the shadow of such objects is curved surface dif- 
fraction. When an incident ray strikes a smooth curved perfectly 
conducting surface at grazing incidence, i.e., at the shadow boundary, 
a part of its energy is diffracted into the shadow region. To 
describe this phenomenon Keller [51] introduced a class of curved 
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surface diffracted rays. These ray paths include the points Q] and 
Q2 which form a curve on the diffracting surface as illustrated in 
Figure 6. The concept of surface rays was actually presented earlier 
in the creeping waves introduced by Franz and Depperman [52,53], 

The results presented in the following discussion are those given 
by Kouyoumjian [54,55]. 



SHADOW 

BOUNDARY 



Figure 6--Diffraction by a smooth curved surface. 


The diffraction by a smooth curved surface is shown in Figure 6 
in which 0 is the source point and P is the observation point. 
Applying Fermat's principle, the Tine 0Q]Q2P is the shortest 
distance between 0 and P which does not penetrate the surface. This 
is a geodesic path with arclength t on the surface. As the field 
propagates along the surface it sheds energy tangentially into the 
shadow region. The preferred coordinate system for this diffraction 
mechanism is composed of $ the unit vector iji the direction of in- 
cidence, £ the normal to the surface S, and ^ the binomial to the 
surface (S = t x n) . The incident field may be resolved into the 
normal and tangential components. It is assumed that these two 
components induce surface ray fields which propagate independently 
of each other. The diffracted field is given by 
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(23) 


/ n \ „ r 1 


E (P) “ E (Q-| ) • [b- ] b 2 u s + sU+s) 


-jks 


where 


r V 2 


s ,h 


N 


dn 


dn 


% ,h (t' )dt' 

_ m 

I e" jkt I D m’ h(Q l )D m’ h(Q 2 )e ^ (24) 


and 


nl,2 n I -j >2 ’ t l ,2 = t h,2* b l ,2 ' b h,2‘ 


The subscripts (superscripts) s,h denote the acoustically soft and 
hard boundary conditions, respectively. The quantity D m (Q) is the 
scalar diffraction coefficient at the point Q and a m is the at- 
tenuation constant which is a function of t because it depends on 
the local radius of curvature and its derivatives. Voltmer [56] 
has obtained diffraction coefficients and attenuation constants to 
order (2/kpg) 2 / 3 (where pg = local radius of curvature of the 
surface along the ray path) which is valid for smooth convex surfaces 
of non-constant curvature. These curved surface diffraction and 
attenuation constants [56] were obtained from an extension of the 
work of Hong [57]. The diffraction coefficients and attenuation 
constants are given in Table I. Note that even though an infinite 
number of modes are indicated in Equation (24), it is generally 
found that sufficiently accurate results for the field in the 
shadow region are obtained by summing on only the first couple of 
modes, provided the field point P is far-removed from the shadow 
boundary. 

The parameters dnl , dn2 and p may be found from differential 
geometry involving the rays and surface; this is discussed at length 
by Levy and Keller [58]. For the important case of cylinders as 
considered in this study, it can be shown that 



and that the caustic distance p is given by 
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Al'(-x) denotes the derivative of A1(-i) with respect to the * 1-90*0 1 
of the Airy function. 




(25b) 


p = s' + t. 


where s' is the distance from the source to the point Q-j . 

The above discussion has been applied to open curved surface so 
far. For a closed surface, each surface ray mode produced at Q-j 
encircles the surface an infinite number of times. In the present 
work, the surfaces are large enough, so that these multiple encircle- 
ments may be neglected. An additional surface ray field must be added, 
however, associated with the shadow boundary at the other end of the 
closed surface to Q^. 

H. Transition Field Associated with Curved Surface Diffraction 


In the analysis of the scattered fields from smooth, convex 
surfaces the pertinent diffraction and reflection coefficients are 
readily available, but they are only valid away from the shadow 
boundaries. At and near the shadow boundaries, the 6TD ray method 
fails. Therefore, in this transition region adjacent to the shadow 
boundary a separate solution is required. 

The general problem of the diffraction of waves by smooth, convex 
surfaces for grazing angles is very complex, and it has been the 
subject of an investigation by various authors. An extensive 
bibliography on various aspects of the problem may be found in the 
papers by Logan [59], and Borovikov and Kinber [60]. The work 
in these references and others are based on the pioneering work of 
Fock [61]. In general, due to the mathematical complexity of the 
solutions, the previous solutions do not appear to be in a tractable 
form for practical engineering problems. 

Pathak [62,63] has modified and extended the transition solution 
in a systematic fashion via ray optical considerations. His result 
is in a form suitable to be used in high frequency scattering by 
arbitrary two-dimensional and three-dimensional smooth convex, 
conducting surfaces. This study will briefly discuss Pathak's [62, 

63] results for diffraction by cylindrical structures. This uniform 
result yields a finite, continuous variation for the total field 
across the shadow boundaries within the transition regions. In 
the lit region, exterior to the transition region the solution reduces 
to the incident and reflected geometrical optics field as presented 
in Section B. In the deep shadow region, exterior to the transition 
region, the solution reduces to the usual GTD surface diffracted ray 
field as given in Section G. The transition problem is illustrated 
in Figure 7a for the lit region (also refer to Figure 2) and in 
Figure 7b for the shadow region (also refer to Figure 6). 
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Figure 7— Curved surface transition function geometry. 
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The transition field for a cylindrical structure can be written 
as 


ft 
E ^ 


E 1 (Q) 


[u l u 2 G s + v l v 2 G h ] 


(26) 


where 

r 


Q = 



for the lit region case 


j^Q-j for the shadow region case 


and G s , Gh are curved surface transition coefficients for the soft 
and hard boundary conditions, respectively . In the lit zone, the 
transition coefficients are given by 



r r 
p l p 2 


-jks 


(p 1 +s)(p 2 +s) 


where 


and 



v 


1 


= e _i 



(27) 


The first term in the square brackets is similar to the Kirchhoff 
edge diffraction coefficient of a half-plane. The function F ( X ) in 
this term is the transition function of the wedge in Section C 
Equation (13). In contrast, the second term in the square brackets 
is dependent on the electrical properties of the surface; this term 
becomes dominant far from the shadow boundary. The functions p*(?) 
and q*U) are universal functions and are the complex conjugates 
of the pU) and q(^) functions defined by Logan [59]. The functions 
e-j ir /4p*(£) and e _ J 7r /^q*(c) are plotted in Figures 8a and 8b, 
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Figure 8a--Plot of e" J 17 ^ p*(0 vs 5 based on Logan's 
tabulated data [59] for pU). 
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Figure 8b. Plot of e -J ^ q*(0 vs K based on Logan's 
tabulated data [59] for q(s). 
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respectively. For a further discussion of the derivation of these 
functions see Reference [59]. The other parameters are defined as 



p (Qn) = longitudinal radius of curvature of the surface at 
9 the point Q 


p, 2 are defined in Section B, Equation (10) for spherical 
wave incidence. 

In the shadow zone, the transition coefficients are given by 



The functions F(x), p*(c) and q*U) and the parameter L are defined 
as in the lit case. The other parameters are defined as. 
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1,2 


k p q^1 ,2^ 


1/3 


? s = 


Q 2 / kD X 1 / 3 

i_fza ] df 


Q-| M 9 


Pn V 2 


p ( t 1 ) is the radius of curvature on the surface at the point 
9 t' on the geodesic path of the ray 


and 


dn-, , dn 0 , p are defined as in Section 6, Equations (25a) and 
1 1 (25b). 
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CHAPTER III 

NEAR FIELD SCATTERING BY A NUMBER OF FLAT PLATES 


A. Introduction 


The Geometrical Theory of Diffraction as discussed in Chapter II 
is applied in this chapter to the near field scattering by a number 
of finite flat plates. This is basically an extension of the 
scattering by a single finite flat plate considered by Burnside [23]. 

In this problem the antenna is in the near field of the three- 
dimensional scattering bodies, and the observation point is in the 
far field. However, it is assumed that any point on the scattering 
body is in the far field of the source. Some assumptions have been 
made about the geometry of the scattering bodies to simplify the 
analysis, that is to minimize multiple scattering mechanisms. 

These assumptions are based on geometries that can be used to obtain 
results for wing mounted antennas. 

B. Near Field Scattering by a Number of Finite Plates 

The near field scattering by a number of finite flat plates 
presented here is a practical application of the Geometrical Theory 
of Diffraction presented earlier. The dimensions of the plates are 
assumed to be large in terms of the wavelength. The geometry of the 
plates is specified by the location of their corners in a general 
coordinate system, as illustrated in Figure 9. The source is 
defined by its location and far-field pattern. The far field 
pattern is appropriate since the source is assumed to be 2 D^/a 
away from the edges of the plates, where D is the maximum dimension 
of the source. The infinitesimal electric and magnetic dipole moments, 
discussed in Section II-D, are the source types considered in this 
analysis. These elements, properly weighted and summed, can be used 
to approximately represent any antenna with a known source distri- 
bution. 

The high frequency scattering by a number of structures in 
proximity to each other is, in general, a very complicated problem. 

All the multiple interactions between the various structures must 
be included to obtain an exact solution to the problem. This of 
course is impractical in most instances. For low frequencies the 
method of moments is used for these types of problems, however, 
for electrically large sized bodies physical limitations on 
computer storage and numerical accuracy make this method im- 
practical. The GTD and its ray optical format is ideal for these 
problems. 
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Figure 9 — Geometry for the near field scattering by 
a number of flat plates. 

In constructing a solution using GTD, one starts by tracing the 
most important rays that will contribute to the total field. It is 
important that all rays be properly shadowed by the scattering 
structures. This, of course, will produce discontinuous fields 
at the shadow boundaries. These discontinuities, however, may be 
compensated, by including the appropriate diffraction mechanisms 
discussed in Chapter II. As the solution progresses, more and more 
higher order terms can be added until the resulting patten is 
relatively smooth. In general , it is impractical to include enough 
higher order terms to produce a completely continuous field. How- 
ever, the discontinuities may be kept to a few dB and in fact the 
general rule is that the actual field normally passes through the 
average of the discontinuity. The resulting radiation patterns are 
usually a good engineering approximation to the true fields. 

The geometry of the multiple flat plates discussed in this report 
is assumed to be such that second order interactions only need to 
be considered. This is a good approximation for the aircraft 
geometries of interest in that the plates are relatively widely 
spaced and do not produce strong corner reflector effects for 
antennas mounted on the wings. The type of interaction terms 
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that have been determined to be important for the geometries in 
this study are listed in Table II. They comprise the major 
contributors to the total scattered field from, a number of flat 
plates. 


TABLE II 

SCATTERED FIELDS INCLUDED IN THE FLAT PLATES SOLUTION 



Symbol 

Description 

Figure 

Equation 

1) 

r 

Simply reflected field 
off a plate 

13 

41 

2) 

E 6 

Diffracted field off 
an edge 

15 

46 

3) 

E 0 

Corner diffracted field off 
a corner 

16 

48 

4) 

r rr 

Double reflected field 
off a plate after re- 
flection by another plate 

21 

49 

5) 

^d 

Diffracted field off an 
edge after reflection by 
a plate 

22 

50 

6) 

r c 

Corner diffracted field 
off a corner after re- 
flection by a plate 

23 

51 

7) 

^dr 

Reflected field off a 
plate after diffraction 
by an edge 

24 

52 


The incident field from the source is given by 

* y\ ^ — j kS 

e\s) = [e s F(0 s ,<j, s ) + ; s G(e s ,* s )] . (29) 


For a far field observation direction given by 

A /N A A 

d = cos<}> s sine s x + sin sin e s y + cos e $ z. 


37 




the phase can be referred to the reference coordinate system by the 

usual far field approximation so that e”^ s /s £ e jkRs*d e'j^d^ 
where d is the distance from the reference point to the far field 
receiver. 

It is essential, in the GTD format, that the source field be set 
to zero whenever the ray from the source to the receiver passes 
through any perfectly conducting obstacle. This can be accomplished 
in many different ways. The scheme used here is to first decide 
whether the ray direction, representing the field., is pointed in 
the same half space in which a plate with normal n is contained. 

This is done by comparing the sign of 

n • (R $ - Ci ) and n • d. (30) 


If they are of the same sign, the ray can not hit the plate. If they 
are of the opposite sign a further test is needed. 

First, the incident ray is projected onto the plane of the plate 
to find its intersection point, as shown in Figure 10, that is 



[n • (R s -C 1 )]d 

A A 

n • d 


(31) 


Now, using an idea based on Cauchy's formula from complex variables, 
that is 



0, no pole in f(z) 
2ttj, one pole in f(z) 


the intersection point can be tested to see whether or not it falls 
within the limits of the plate. This is illustrated in Figure 11. 
It is easy to show that 


e 


m 


tan 


(l V R t> * (WV 

<w • <w*t> 


(32a) 


which leads to the test, if 
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X 


Figure 10--Iritersection of scatter direction ray with plate. 
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(a) RAY HITS PLATE 


m + I 



(b) RAY DOES NO' 


<tt , no hit occurs 
>tt , a hit occurs 


(32b) 


M 

I 


m=l 


m 


All the plates may be tested in this way to see if they shadow 
the ray. 

The singly reflected field from a flat plate is considered 
next. The first step in this calculation is to find the location 
of the image source, which is uniquely determined once the plane 
of the flat plate is defined by its normal. In fact, the image 
is located along a line which is orthogonal to the plate and 
positioned an equal distance on the opposite side of the plate. 

This location can be found analytically using the geometry illustrated 
in Figure 12. The image position vector (Rj) can be determined using 
the following expression 

Ri = R s - 2[n • (R s - )]n. (33) 


With the image position known, one needs to determine if the 
reflected field contributes to the total scattered field. The 
reflected field exists if the ray from the image source in the scatter 
direction passes through the finite limits of the plate. The same 
£rocedure outlined in Equations (30-32) is used for this except 
Ri replaces R s . In this case if a hit on the plate does occur, 
reflections could exist. One must also test that the incident 
ray and the scatter direction ray are not shadowed by another plate. 
The incident ray (I) may be found by reflecting the scatter direction 
ray into the plane of the plate. This is given by the expression 

I = d - 2(n-d)n . 

The incident angles for the source may be found from 



where 

' = V + v + 
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Figure 12— Geometry for location of image source. 


or 


I = cos <f>.j sin x + sin <}>.. sin y + cos e ^z. 


Once it is determined that the reflections do occur, the value of 
the reflected field can be found. For the jth plate the reflected 
field from the image source can be written as 


E J(W = C 0 S F r (V* s ) + ^ S G r ( e s ^ s )] - 


•jks 


(34) 


using the geometry illustrated in Figure 13. The vectors 0 s and <l> s 
are related to the scatter direction d and are given by 


0 s = cos<t> cos© x + sin<S> cos0 v - sin0 z 
s s s s s 

A A A. 

4> s = -sin<i> s x + cos<l> s y. 

The incident field is given by 

a # „ j 

E^s 1 ) = [0 i F(0 i ,4» i ) + ♦ 1 6(0 1 ,* i )] 

where 

/N . A A A 

0 1 = cos^. cosQ^x + sin'i'^cosQ^y - sin0..z 


(35) 


(f) 1 = -sin^^x + cos^y . 

A tangent vector (t) can be found on the plate such that 
2 _ n x d 

L “ A A • 

n x d 


(36) 


The reflected field can now be found from the boundary conditions 
that must be satisfied on the perfectly conducting flat plate 
surfaces. They are given by 


E (on plate) = n 


E 1 


(on plate) 


(37a) 
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X 


Figure 13 — Reflected field geometry. 


44 


(37b) 


E r 


(on plate) = - t 


E 1 


(on plate) 


In matrix notation Equations (37a and 37b) can be written, using 
Equations (34 and 35) and the fact that s=s' on the plate, as 



where 


! 

< -e- 

< c 

•i— 

< CD 

< C 
l 


r~ 

U_ 


1 

c 

< 

1 

A A J A A a 

_t • 8 1 t • 


CD 
J 

II 

\ 

4-> 

c 

1 


(38a) 


(38b) 


Now solving the above expression for the unknown reflected field 
pattern functions gives 


f r = (t-j S )A n + (n-j^A 1 

* A A »» A A p A A _ A A — 

(n.e s )(t.^ s ) - (n-(j> s )(t-9 S ) 


(39a) 


G r = 


■[(n e s )A t + (t.e s )A n ] 

A A m A A * A A A A A a ^ * 

(n-e s )(t- 4 » s ) - (n.<}, s )(t. 0 s ) 


(39b) 


Thus, in the far field the reflected field from the jth plate can be 
written using Equations (34 and 39) as 


■ft 


^(e s ,<j, s ) = [e s F r (e s ,<j> s ) + ^G r (e s ,<|> s )]e 


:s 


jkR i -d e -jkd 


. (40) 


The total reflected field from the P flat plates is, then, given 
by 



remembering that if the incident or reflected field is shadowed 
by any pi ate,, the reflected field for that plate is taken to be 
zero. Also, if the source is directly mounted on the surface of a 
plate, such as a slot antenna, there is assumed to be no re- 
flections from that plate. 

The shadowing of the incident and reflected rays by the flat 
plates is compensated for in the total field by the addition of the 
three-dimensional wedge diffracted field of Equation (16). It is 
known that for a given scatter direction there is only one point 
along an infinitely long straight edge at which the diffracted field 
can emanate for a near zone source. Thus, one needs to determine 
whether a diffraction point lies on the finite length of each of the 
M edges on the P plates so that it can be included in the total 
scattered field. Note that it is possible for two plates to have 
a common edge. The common edges are treated as one edge with a wedge 
angle other than n=2 and therefore it must only be included once in 
the scattered field. 

To determine if a given edge has a diffraction point that will 
contribute to the total field in the scatter direction, bounds on 
the possible diffraction cone angles g 0 can be obtained for each 
of the edges forming the P plates. Diffractions will occur if the 
following inequality is met 


(C -R ) 
m s' 

C-Rs 

m 


~ Z ~ < C m + r R s> 

e < d ■ e_ < . e 


m 


m 


C m+1 " R s 


m 


(42) 


otherwise, diffractions for that edge will not occur in the given 
scatter direction. The diffraction point (Rd) for an infinite 
edge and an observation point, as illustrated in Figure 14, can 
be obtained in a straight-forward manner. For a given edge, a 
perpendicular vector from the source to the edge can be found 
from 


t = (R - C ) • e , 
p s m m 


R = C + te . 
p m pm 


(43a) 

(43b) 


The scatter direction is known so that 
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(43c) 


d • e 


COt 6 = 


m 


O I ~ I o 

V l-(d-e ) 2 
m 


The diffraction point can now be found from 


Rj = R + s„ cot e m 
d p p o m 


where 


s = |R - R | 
p's p 1 


(43d) 

(43e) 


and the incident vector is given by 



As before in the discussion of the reflected field the incident 
field and diffracted field must be checked to see if they are shadowed 
by any other structure. This can be done in the same manner given in 
Equations (30-32). If either field is shadowed the result is set 
equal to zero for that edge. 

Once the diffraction point is known, the diffracted field can be 
found from Equation (16). Since the source produces a spherical wave 
incidence and the scatter direction is in the far field (S-h*>), the 
diffracted field (with its phase referred to the reference coordinate 
system) is given by 

eJ (Q e ) 

e2 (q e ) 

jk e -jkd 

e ~~T~ (44) 



— 1 
m 

: CL 


-D - ir- 

S 9 

0 

11 


s jk 

9(f>' an 

r 

m 

CL 

1 

u 

0 


.0 .L5.L 

h jk 9 <j>' an J 


where 


L = s 


. 2 . 
sin e r 


(45a) 
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r 


n w = wedge angle number 


1 tan' -tt' 


^n. -n •> 

• J 


for half plane edges 

for angle between ni 
and n. plate with 

SJ 

plates i and j 
having a common edge 


and 


(45b) 


e!(Q e ) = E i (e i ,4> i ) • (V 


(45c) 


e](Q e ) = E 1 (e. ,<J> i ) • 4)' 


s' = I 


(45d) 

(45e) 


0 i = tan 


-1 'J 




(45f) 


*1 


tan 


-i r* 


(45g) 


The incident field is given by 

i /s ^ — j ks 1 

E ( 9 ^ ) = [0 n 'F(0. ,4>.) + ,<|>.)] “n — • 

In order to define the diffraction angles it is necessary to 
define a coordinate system at the diffraction point a§ illustrated 
in Figure 15. The three orthogonal unit vectors are e m , the unit 
vector along the edge; n, the normal to the surface; and em=nxe m , 
the vector perpendicular to the edge and lying in the plane of the 
flat plate. The incident vector direction can now be written in 
terms of this coordinate system as 
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Figure 15--Edge coordinate system at point of diffraction. 
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Using the above expression one finds that 



The scatter direction unit vector can, also, be decomposed in 
terms of this coordinate system as 

d = e ' cos<p sin6 + n si n<J> sinB + e„ cosB„ 
m o o m o 


which results in 


<{> = tan 


-1 1 d • n 


The half-cone angle of the diffracted cone of rays is found from 

1 A A 

= sin - Id x el 
o 1 m' 


and the incident cone angle is found from the laws of diffraction to 
be 


The vector directions used to determine the diffracted field are 
defined by the following expressions: 

<)>' = -e^ sin<j>' + n cos<f>' 

A A A 

<f> = -e^ sin<|> + n cos <}> 

AAA A A A 

B' = <J>' x I and B„ = <|> x d. 
o o 

The total diffracted field from a general mth edge of the jth plate 
can now be found to be 



E*? (9 ,<j> ) = E 1 

j,m v s’ T s' 


*5 0 + E d i 


The total diffracted field is the superposition of all the diffracted 
fields from all the edges on the P plates that exist in that direction 
and is given by 







(46) 


The diffraction from an edge as discussed above is based on an 
infinite edge. The plates of interest here have finite length edges. 
The fact that the straight edge diffraction point is lost past the 
corners of the edge produces a discontinuity that must be compensated 
for in some way. In the Geometrical Theory of Diffraction this is 
accomplished using a corner diffraction term. This diffraction 
mechanism is very small in magnitude except when the source is close 
to the corner or if the radiation pattern is taken directly across 
the corner. For wing mounted antennas, this situation occurs fre- 
quently enough to warrant the inclusion of this term. As discussed 
in Chapter II, Section F, the corner diffraction term is available 
at present as a heuristical ly derived function. This function has 
proven to be useful as will be seen below. 

The corner diffraction term is assumed to be localized to a 
small region around the truncation point of the edge, that is, the 
corner. For each edge on the plates, there are two terms associated 
with each truncation of the edge. This is illustrated in Figure 16. 
The diffraction point is found for a given scatter direction using 
Equations (43a-e). This time, however, the diffraction point is 
not bound by the finite limits of the edge. It may fall off the 
edge on an imaginary extension of the real finite edge as shown 
in Figure 16. Since this term will be small away from the edge, 
the limits for including the terms associated with the two trun- 
cations are given as 


< C m- R s> 



e m - 0.1 < d 
m 



(C ) 
v m s ' 

C m -R_ 
m s 


v 0 - 1 • 


for the first (mth) corner 
of the mth edge. 


and 
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< W R s> 
SlH-1 " R S 


'm 


0.1 < d • e < . e +0.1 

“ ID 


m 


C m+r R s 


for the second (m+ist) corner of 
the mth edge. 


This bound has been set heuristical ly , but it appears to give good 
results in the cases tried. Once the diffraction point is found, 
tests are made on the incident direction and scatter direction 
to see if they are shadowed by any obstacles as was done before. 

After it is determined whether the corner terms will contribute 
to the scatter direction, the terms for the first and second edge 
truncations of the mth edge are computed from Equation (22a), that 
is 



. ir 


_J 4 

sin£ c e 

2tt(cos8 -cosB 
' 0 c 


7 


F [ks c a(6 0 -B c )] 


jk R d -d 


e "jkd 

d 


(47) 


for a scatter direction in the far field. The quantities I and M 
are defined in Equations (22b and c). 


For the jth plate, the mth edge, and the ith truncation (i=l or 2) 
the corner diffracted field can be written as 




J ,m,i 


i( 0 e-i) = E l 


B o + 


eS 


0. 


All the variables are defined as in the discussion for the edge dif- 
fraction term defined previously, except for the terms that follow. 
For the first truncation point (mth corner), 
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COS g_ = 


C |R_ 


m 


z‘ = C m - R. 
m d 


For the second truncation point (m+lst corner), 


R c = R s - C m + 1 


s c = l R c' 


cos 6„ = — — • (-ej 


C !r_ 


m' 


2 ‘ = C m+1 " R d 


The corner diffraction term is not presently applicable to corners 
that are formed by the junction of two plates. The corner term is 
neglected when this case arises. The total corner diffracted fields 
of all the corner diffracted fields that exist in the scatter 
direction is given by 



p "j 

i v 

j=lm=l 






(48) 


The four terms discussed above comprise all the mechanisms 
needed to determine the scattered fields for a single plate. The 
total radiated field in the presence of a single plate is given, 
using Equations (29), (41), (46), and (48), by 

E ( 9 $ ) = E i (e $ ,(i> s ) + ( 0 s ,^ s ) + E^(0 s ,<|> s ) + ( © s ,4> s ) • 
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This solution has been compared with measured results using a short 
dipole antenna mounted above a flat square plate [23]. This is il- 
lustrated in Figure 17 for a dipole mounted above the center of the 
square plate. To illustrate the effect of including the slope dif- 
fraction term, Figure 18 compares the radiated field of a short slot 
mounted on a plate with and without the slope diffraction term. 

The effects of the corner diffraction term is shown by obtaining 
the radiation pattern across the corners of a square plate using a 
short monopole mounted at the center of the plate, as shown in the 
insert of Figure 19. The radiation pattern with and without 
corner diffraction is compared in Figure 19. Note that the 
discontinuity is completely eliminated. The radiation pattern 
with corner diffraction is compared with an experimental result in 
Figure 20. The good agreement in the comparison would tend to confirm 
the approximation made in the corner diffraction coefficient. 

If more than one plate is desired, the interaction between 
the plates must be included. These higher order terms are obtained 
in a step by step manner similar to the analysis carried out for 
the mechanisms discussed above. First, the double reflected fields 
will be discussed. The double reflected field, as illustrated in 
Figure 21, can be viewed as a field emanating from a second image 
obtained by reflecting the first image of the jth plate into the 
kth plate, that is 

*11 = *1 ' 2 tn k * ( *i " ^kl^ n k 


where R-j is obtained from Equation (33) for the jth plate and nj^ 
refers to the normal of the kth plate. The ray emanating from Rjj 
in the scatter direction 3 must be tested to see if it hits the kth 
plate. This is done using the procedure outlined in Equations 
(30-32). If a hit does occur the analysis proceeds to the next 
step which is to find the incident ray 

I k = d - 2(n R • d)n k . 

The ray emanating from R-j in the direction 1^ must be tested to see 
if it hits the jth plate as was done above. The incident field on 
the jth plate is then found from 


Of course, each of the rays must be tested to see that they are 
not shadowed by any other plate. 
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Figure 17--Ee radiation pattern for a small dipole mounted 
above a square plate for e s = 90° and 0£<j> s <360° 
at freq. = 10.43 GHz. 
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WITH SLOPE 
DIFFRACTION 



Figure 18--Comparison of the E 0 radiation pattern of a strip slot 
on a plate with and without slope diffraction. 
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Figure 19--Comparison of the ; E e radiation pattern of a short 
monopole on a square plate (A=4x) with or with- 
out corner diffraction taken in the plane $ =45°. 
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CALCULATED 

MEASURED 



Figure 20--Comparison of the calculated and measured Eg rad 
ation pattern of a short monopole on a square 
plate ( A=5 . 51a) taken in the plane <f>=45°. 


60 




The doubly reflected field can now be obtained, if all of the 
tests determine that a field exists. First, the reflected field is 
obtained from the jth plate as was done in Equations (34^40), 
except with e^. . substituted for © ,4> where 


I 


k 


J kx* + V 


+ I 





This means that the reflected field from the jth plate, which is 
the incident field on the kth plate, is given by 



ki 


= [ 0 k Fr(e ki ^ 


;’G r (e kf 




e-J k5 ' 
s ' 


The boundary conditions on the kth plate are satisfied such that the 
doubly reflected field is found. In matrix notation this can be 
written as 


'n k -0 S 

n 2 s 1 

"k • + 


■ F rr- 


A n 

/V A _ 

Lv 9 

A A - 

V* - 


-G rr . 


--A 1 . 


where 
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"v®k 

Y*k' 


-F r - 


' A n * 

L t. - e? 

A A * 

V*i, - 


. G r _ 


. A* . 


Thus, solving the above expression for the unknown reflected field 
gain functions gives 

(t k .; S )A n + (n k 4 S )A t 
F " (n k -i s )(t k -i s ) - <n k 4 s )(t k -e s ) 

^ rr _ -[(n k -e s )A t + ( t k -e s )A n ] 

(n k -0 S )(t k -i S ) - (n k 4 s )(t k -e s ) 


In the far field the reflected field from the kth plate after its 
original reflection from the jth plate can be written as 


^(e s ,+ s ) = [eV r (e c ,»J * ;V r (e.,*J]e jkK ir a £ 


-jkd 


s’ T s' 


s' T s' 


The total double reflected field from the P plates is, then, given 
by 

* I, j (49 > 

j^k 

— rr • 

for all the Ej k fields that are not shadowed. Note that there is 

no double reflection off of the same plate, i.e., when j=k. 
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The fact that the kth plate has finite limits is compensated 
by a diffraction term off the M|< edges of the kth plate. Thus, the 
reflected field from the jth plate must create diffracted fields from 
the edges composing the kth plate. This reflected-diffracted term 
is illustrated in Figure 22. In this case, as in the_previous case, 
the reflected field from the jth plate and its image Rj may be sub- 
stituted into the needed equations for the incident field and 
source point, respectively. That is, bounds for the diffraction 
point on the mth edge of the kth plate are obtained from 




: j : ^ c k,m+i R P ^ 

e < d • e < . e 

m m itt k i m 


|C 


k,m+l ~ R i 


as in Equation_(42) . The diffraction point Rd is found from Equation 
(43a-e) , with R-j substituted! for R s . The incident vector on the 
edge is given b Y ^ = R d " R i' 
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Once the existence tests on the rays are performed, the re- 
flected-diffracted field can be found from 


m 

Q- 

1 


1 

o 

in 

0 
i 

1 


_E X>' 






-tf- 


Lo -D h - 


L m 

1 5 

O 

m 


s' 

\ 


Jk Rd • d 



Note that the slope term has not been used here. This is due to 
the fact that the slope term did not appear to be a significant 
contributor to the fields for the geometries of interest. The 
incident field (reflected field from the jth plate) on the edge 
is given by, Equation (40), 


E^(Q e ) K^ld^ki) + tc^ki »+ki^ 


e _jks ' 

~s r ~ 


in which case 


e m(Q e ) E j( 0 k1’*ki) * 


E I (Q E } = E j (0 ki’ <t> ki ) • ♦' 


s' = 


It 


e. • 

ki 


= tan 


2 2 

I + I 
*kx K.y 


I 


kz 


* ki = tan ’fe 


All the other parameters are as determined previously. 
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The reflected-diffracted field from the mth edge of the kth 
plate is then given as 


|rd 


j»k,m 




The total reflected-diffracted field of all the reflected-diffracted 
fields that exist is given by 


F d (e s ,4 s ) 


k -F-rd . x 
III E j»k ,m' e s ’^s 
j=l k=l m=l 
j/k 


(50) 


The same procedure for the reflected-diffracted field is applied 
to the reflected-corner diffracted field so it will not be discussed 
in detail here. This field is illustrated in Figure 23. The 
corner diffracted term for the mth edge on the kth plate is found 
from Equation (47) with the reflected field of the jth plate given 
by Equation (34) used as the incident field. The total corner dif- 
fracted field in the scatter direction is given by 


r c ( 0 c ,* ) = 


P P M k 2 

l l l l 

0=1 k=l m=l i=l 
j7k 


E j!k,m,i (0 s 5<f, s } 


(51) 


The diffracted field of the mth edge on the jth plate can be 
reflected by another plate (call it the kth plate). This diffracted- 
reflected field is, in general, of the same order of magnitude as 
the reflected-diffracted field discussed above. This field type 
is illustrated in Figure 24. 

The diffracted point for the mth edge on the jth plate, if it 
exists, is foynd using Equation (43a-e), where must be sub- 
stituted for d. The diffraction direction in this case is the 
incident direction for reflection of the kth plate and is given 
by 


A A A A A 

I k = d - 2(n k -d)n k . 


Once the diffraction point R . is found the imaged diffraction point 
is given by 



= # d - 2 <V R d>V 
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it* CORNER 



Figure 23 — Reflected-corner diffracted field geometry. 
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The imaged diffraction point is then used to see if a reflected 
field exists by virtue of an incident ray hit on the plate. The 
diffracted field can then be found from Equations (44 and 45a-g) , 
with the appropriate incident direction for the incident field and 
diffraction direction for the diffracted field. In this case, the 
slope term is not included. This diffracted field can now be used 
in Equations (38-40) as the incident field for the reflection 
mechanism. The details will not be given here. 


In the far-field the reflected field of the kth plate from 
the diffracted field of mth edge on the jth plate is given by 


= [9 S F dr (e s .4» s ) ♦ ; S G dr (6 s ,* s )] ( 


jkR 


di d e' jkd 


The total diffracted-reflected field of all the diffracted-reflected 
fields that exist is given by 




P M. P 

I I I 

j=l m=l k=l 
jVk 


E^ r . (e ,<t> ) 

j,m,k v s ,M V 


(52) 


where Ej „ k = 0 if the incident, diffracted, or reflected ray is 
shadoweab^ any obstacle. 

In certain cases, another higher order term can be significant. 
This is a double diffracted term, illustrated in Figure 25. This 
term is normally small unless the two edges and the source are 
close together or one of the edges lies in the shadow boundary of 
the other edge. In general, this is a fairly complicated ray tracing 
problem without producing results that are significant to the total 
field. When the two edges lie in the same plane, the problem 
simplifies omewhat; however, the results are still not significant 
except for small angular regions about the shadow boundary. The 
problem as to when it is necessary to compute this term for a 
general structure, is under investigation. The double diffracted 
term has not been added at this time and there has not been sub- 
stantial deterioration of the patterns. 
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Figure 25--Double diffracted field geometry. 


The total field for a source in the presence of a number of 
plates can now be approximated by superimposing the individual 
scattered fields listed in Table II, as 

E(e s „g - r + r + f + r + r r + r d + r c + r lr . m 


This represents the sum of the most significant terms for a number 
of plates that are fairly widely separated. In the case of plates 
with common edges, the angle between the plates should be great 
enough not to produce a large number of multiple reflections. 

There is the possibility that other terms will be significant for 
particular geometrical configurations. These are seen as dis- 
continuities in the pattern. However, using the above prescription, 
these rays may be identified and included in the total field. An 
engineering judgment must be made as to whether the time expended 
in obtaining these higher order terms will produce a corresponding 
benefit in the final result. 
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This solution can be used in a wide variety of ways. To il- 
lustrate a simple case, this solution is compared with results 
obtained using the newly developed moment method patch technique [7] 
as illustrated in Figures 26a and b and 27a and b. Note that for 
these comparisons the plate dimensions are quite small in terms of 
the wavelength to accommodate an overlap between the low frequency 
patch solution and the higher frequency GTD solution. This accounts 
for the discontinuities in the patterns using the present GTD 
solution. 

More complicated structures can also be analyzed using this 
solution. For example, the scattering from a building or a ship 
can be modelled using flat plates. In this study, the plates 
are used as models for the wings and horizontal and vertical 
stabilizers of aircraft. They could also be used to model the 
wing flap, or the stores and engines on the aircraft. Even more 
complicated structures such as the fuselage of the aircraft can 
be modelled in a crude way. Such a model was constructed for an 
F-4 aircraft with a slot mounted on its wing before the more 
advanced model using an elliptic cylinder discussed in Chapter V 
was developed [22]. 

The fuselage of an F-4 aircraft is modelled as a lopsided flat 
plate as illustrated in Figure 28. The wing is modelled as in 
Figure 29, with a 17.5° bend between the two sections about the 
common edges numbered IT and 13'. In order to check the validity 
of the approximation made in the flat plate model of the F-4, 
experimental measured patterns were taken at the Naval Weapons 
Center using a scale model of an F-4 aircraft, shown in Figure 30. 
Four different pattern cuts were made as illustrated in Figures 
31a-d. Five different antenna locations on the wing were used 
as shown in Figure 32. The antenna was a single cavity-backed 
slot of length 0.748A in free space with only the first waveguide 
mode excited. The comparison of experimental and calculated main 
beam elevation plane pattern results (refer to Figure 31a) for 
the slot placed at a 45° angle with respect to the fuselage axis 
at locations #2 and #4 are shown in Figures 33a and b, respectively. 
The pattern in the main beam roll plane (Figure 31c) for an antenna 
at location #2 is shown in Figure 33c. The elevation plane pattern 
(Figure 31b) for slot antenna location #2 is shown in Figure 33d 
and its roll plane pattern (Figure 31a) is shown in Figure 33e. 

At the time of the comparisons the slope diffraction, corner dif- 
fraction, and double diffraction fields had not been added. The 
resulting discontinuities are marked on the patterns for illustrative 
purposes. The overall character of the measured and calculated 
patterns appear to agree; however, the finer details of the patterns 
do not compare very wel 1 . 
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Figure 26a--Comparison of measured, reaction and GTD calculated 
H-plane pattern results for a x/Z dipole near a 
corner reflector. 
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Figure ,26b--Comparison of reaction and GTD E-plane pattern results 
(ct=90°) . 
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Figure 27a--Comparison of reaction and GTD H-plane pattern results 
(a=270°) . 
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Figure 28--Fuselage geometry for best plate model 
of F-4 aircraft. 


75 






















10 


15 


Z 



POINT 

z 

KU 

10 

-1.840 

E EM 

1 1 

-1.840 

0.427 

12 

-0.240 

1.868 

13 

0,187 

1.868 

14 

-0.050 

0.369 

15 


EU 

■DU 

-0.876 

1.295 

m 

0.0966 

1.295 


Figure 29--Wing geometry of F-4 aircraft bent plate model. 


6 





















77 


Figure 30--Naval Weapon Center's F-4 aircraft model used in measurements. 






Figure 32--Antenna locations. on F-4 wing used in 
measured results. 
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Figure 33a--Main beam elevation plane pattern of single slot 
mount on F-4 model at location #2 compared 
with bent plate aircraft result. 
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Figure 33b--Main beam elevation plane pattern of single slot 
mounted on F-4 model at location #4 compared with 
bent plate aircraft result. 
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Figure 33c--Main beam roll plane pattern of single slot mounted 
on F-4 model at location #2 compared with bent plate 
aircraft result. 
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CALCULATED 



Figure 33d--Elevation plane pattern of single slot mounted on 
F-4 model at location #2 compared with bent plate 
aircraft result. 
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Figure 33e--Roll plane pattern of single slot mounted on F-4 
model at location #2 compared with bent plate 
aircraft result f 
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Figure 33f — Elevation plane pattern of single slot mounted on 
F-4 wing facing aft compared with bent plate 
aircraft result. 
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The above comparisons, especially the roll plane result in 
Figure 33e in which scattering from the fuselage plays an important 
part, indicate that the flat plate representation does not adequately 
represent the smooth curved surface of the aircraft fuselage. This 
leads to the development of a more complex fuselage representation 
presented in the following Chapters. To minimize the fuselage 
effects on the flat plate fuselage, a measurement was made on the 
F-4 model with the antenna main beam pointed straight behind the 
aircraft. This result is shown in Figure 33f. Some of the detail 
compares slightly better; however, the shifting of some of the 
ripple and the broadnesss of the pattern at 6 = 180° indicates a 
problem with the measurements that turns up in all the comparisons. 
The problem is that the measurements were taken in the near field 
of the aircraft. As a result, one cannot directly compare the 
measured near-field and calculated far-field patterns. Of course, 
there is also the possibility of slight inaccuracies in the measured 
results due to variations in placement of the antennas, ground 
reflections, etc. The comparisons are good enough, however, to 
show that the flat plate model is moderately successful although 
it needs improvement for aircraft fuselage modelling. 

The question as to whether a flat plate adequately models a 
wing with finite thickness and a rounded edge with an antenna 
mounted on its surface is studied [32] by comparing the results of 
a two dimensional thick wing model with the knife edge used in 
the flat plate model. The rounded wing tip is modelled by a 
circular cylinder mounted on a thick wall as shown in the insert 
of Figure 34. The radiation patterns are compared in Figure 34 for 
a tip of 1/16A radius and the knife edge model. Note there is 
very little difference in the two results. This implies that the 
knife edge is a good approximation for a wing tip at the frequencies 
of interest. 
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Figure 34 — Comparison of the radiation patterns of a slot on 
a thick edge plate to that of a slot on a knife 
edge plate. 


87 



CHAPTER IV 

NEAR FIELD SCATTERING BY A FINITE ELLIPTIC CYLINDER 


A. Introduction 


The near field scattering by a finite length, perfectly- 
conducting elliptic cylinder is analyzed using the Geometrical 
Theory of Diffraction as discussed in Chapter II. The near field 
scattering from a finite circular cylinder was studied by Burnside 
[23] using a combination of modal solutions, physical optics, and 
equivalent currents. The present work is now possible due to the 
curved wedge diffraction coefficients of Kouyoumjian and Pathak 
[34] and the curved surface transition function solution of Pathak 
[62,63]. 

In this study, the antenna is in the near field of the three 
dimensional scattering body, but it may not be on the surface of 
the cylinder. The scattering body is assumed to be in the far 
field of the source. The observation point is also in the far 
field. 

The finite elliptic cylinder is of interest because it can be 
used to simulate aircraft fuselages, engines, and stores. The 
elliptic cylinder geometry allows for variations in its dimensions 
such that a best fit can be obtained for a variety of actual air- 
craft fuselage shapes. 

B . Near Field Scattering by a Finite Elliptic Cyli nder 

The finite length, perfectly conducting elliptic cylinder 
analyzed in this study is composed of a cylinder with elliptic 
cross-section truncated by two flat end caps, whose normals lie 
in the x-z plane such that they can make an arbitrary angle with 
the cylinder axis. The end cap cylinder junction forms a curved 
wedge. The dimensions of the structure are assumed to be large in 
terms of the wavelength. The geometry is illustrated in Figure 35. 
The source is defined by its location and far-field pattern. The 
source is assumed to be an infinitesimal electric or magnetic dipole 
moment. 

The incident field is given in Equation (29), and as before the 
scatter direction is given as 


A A A /> 

d = cos <f> sin e s x + sin 4 > s sin e s y + cose s z. 


(54) 
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As mentioned previously, it is important that all fields be set to 
zero if they are shadowed by any obstacle. In this case, the cylinder 
is assumed to be the only structure present. 



Figure 36--Geometry for tangent directions to elliptic cylinder. 

The method used to determine if a ray is shadowed or not by the 
cylinder is based on limits defined by vectors emanating from the 
source such that they are tangent to the cylinder surface as shown 
in Figure 36. These vectors til and ti 2 define the top and bottom 
shadow boundaries of the source, respectively. The coordinates of 
each tangent point associated with the source shadow boundaries can 
be found by setting an incident vector from the source equal to the 
general unit tangent vector to the elliptic surface. This is 
written as 
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( 55 ) 


I = +t 
xy e 


• A 
where the minus sign is associated with the top shadov( boundary (tj*]), 

and the plus sign is for the bottom shadow boundary (t-j 2) • The 

general unit tangent vector to the elliptic surface at the point 


— ^ A 

R = a cos vx + bsinvy+zz 
e J e 

on the surface is given by 



As a unit tangent vector this is given by 




a 

sin v x + b cos v y 

2 2 2 2 
sin v + b cos v 


(56a) 


The tangent vector need only be defined in the x-y plane, so the 
vector emanating from the source is defined only in the x-y plane 
here as 

A A 

(a cos v - x)x + (b sin v - y )y 
I xy - = — ... <56b) 

Ay / 2 2 

\/(a cos v - x $ ) + (b sin v - y $ ) 

where the source point is given by 

R" = x 0 x + yv + z„z. 
s s s 

With the substitutions 

x t y t 

cos v = — and sin v = -e 
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in Equations (56a) and (56b) to simply the notation ,, Equation (55) 
can be used to solve for the tangent points (x f , y., z f ) when it 
can then be used to solve for the tangent points (x. , y., z.) when 
it is written as t t t 


(V x s> 




W (y t -y s ) 2 J^yfbV t 


(57a) 


and 


7(x t -x s ) 2 +(y t -y s ) 2 


+b 2 x. 


J7&v t 


(57b) 


This is accomplished in the following manner. First, dividing 
Equations (57a) by (57b) and rearranging the result gives 



Vs . Vs 
— — - + — — . 


Then, noting that 



(58a) 


is the equation of an ellipse, this leads to the fact that 


Vs 


y t y s 

+ -V = 1 


(58b) 


The y-coordinate of the tangent point can be written in terms of 
the x-coordinate from Equation (58b) as 
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1 


(58c) 




Then, substituting Equation (58c) into Equation (58a) yields the 
quadratic equation on the x-coordinate 



Hence, this can be solved for the x-coordinate of the tangent points 
so that 


2 2 2 7 ~ 2 2 2 2 ? 

1,2 

* <a 2 y 2 ; b 2 x 2 ) 


(59a) 


The y-coordinate of the tangent points are found by substituting 
this into Equation (58c) which gives 


A 2 

1,2 _ d b J s 
y t 


a"b"y c - b 2 xja 2 y 2 +b 2 x 2 -a 2 b 2 


(a 2 y 2 + b 2 x 2 ) 


(59b) 


The top tangent point (#1) is defined with the minus sign in 
Equation (59a) and the plus sign in Equation (59b). The bottom 
tangent point (#2) is defined with the plus sign in Equation 
(59a) and the negative sign in Equation (59b). Then 


v il ,2 tan 



(60) 


so that the tangent vectors are 


t.-, = - 1 
il y e , -v=v.i 


t. 0 = t I 
i2 e'v = v *2 

A 

where t is defined in Equation (56a), 
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The shadowing algorithm follows such that if 


t 


il 




or 




< 



( 61 ) 


where d xy = cos<j> s x + sin 4> s y 

the ray in the d X y direction does not hit the cylinder, that is 
shadowing does not occur. This is illustrated in Figure 37. 

If Equation (61) is not satisfied a further test must be made due 
to the finite length of the cylinder. The incident vector i X y, 
Equation (56b), is assumed to be between the t-j -| and t-j2 directions. 
The parameter v at which the incident vector would hit the cylinder 
in two dimensions is found from 


This leads to the equation 


sin v 


-a sin<f> s cos v + (x g sin <j> s - y g cos<f> s ) 
~ -b cos<(> s 


(62) 


2 

Squaring both sides of this equation, and substituting in sin v = 
1 - cos^v and rearranging gives the quadratic equation 


2 2 2 2 2 

(a sin <(> + b cos <j> )cos v - 2a si n^ (x sin -y <j> ) cos v 

O J J O J J J 


2 2 2 
-b cos <|> s + (x s sint)) s -y s cos^ s ) = 0 


This may be solved for cos v, then cos v may be substituted into 
Equation (62) to find sin v. The sign of the quadratic solution 
is picked such that the |I X y| is a minimum. The z-component of the 
three dimensional incident vector 


A A A 

I = (a cos v - x $ )x + (b sin v - y $ )y + (z e -z $ )z 


(63) 


can be found from 
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cos -1 (t i 2 *d xy ) 





cose„ 

Z = Z + 7 

e s cos<}> sine. 


(x e - x s ). 


where x e = a cos v. The incident ray hit the cylinder and is 
shadowed, as in Figure 38, if 

z” + x„ cote" < < z* + x cot e + , 

c e c e c e c 


where 


z" = origin of 


posi tive 
negative 


end cap 




= angle of 


A 


positive 

negative 


end cap with z-axis, 



Figure 38- -Geometry of intersection point of ray 
in three dimensions. 


96 



Once it is determined the incident field is not shadowed, the far 
field incident field is given by Equation (29). This shadowing 
algorithm can be used for any point of origin and type of 
interaction ray that follows in Chapter V. 



Figure 39 — Geometry for determination of reflection point on 
elliptic cylinder. 

In order to find the reflected field off A the curved surface 
of the cylinder in a given scatter direction d, as shown in Figure 
39, the reflection point must first be found. ^ This is assuming 
that 3 is not in the shadow ^region, that is, d does satisfy the 
test in Equation (61). If d is in the shadow region the reflected 
field does not exist. To find the reflection point in the lit 
zone, the laws of reflection as given by Equations (4a and b) are 
the starting point. These equations are now written such that 

n . I _ |?| _ n x T 

n .. cT }TJ n" x d” 
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or 


(n x I)(n-d) + (n-I)(nxd) = 0. 


(64) 


Note that the normalization factors may be removed. This will 
simplify the solutions to follow. The scatter direction given 
in Equation (54), and the incident vector given in Equation (63) 
may be written as 


■ d x ( V*s )x + yv+s )y + d z (6 s> 2 


and 


I ( v ,z ) = I x (v)x + I y (v)y + I z ( z ) z . 


The normal to the curved surface is given in general by 


n = t x z = b cos v x + a sin v y 


(65a) 


or 


A A 

n(v) = n x (v)x + n (v)y. (65b) 

Performing the necessary dot and cross products, Equation (64) 
reduces to two equations. 


<Vx + "y , y )(n x d y'"y d x ) + <"x d x + n y d y )(n x I y- n y I x ) = 0 (66a) 
and 

( "x‘x + Vy> d z + ("x d x + "yV'z = °- < 66b > 

For a given scatter direction, the location of rT is usually not 
known a priori. It has to be solved for using the above two equa- 
tions. Actually, Equation (66a) contains only one unknown, that is 
the parameter v. Once this is found the appropriate v may be 
substituted into Equation (66b) to find the unknown z . 
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Let Equation (66a) be written as 
3 1 B 2 + 6 3 e 4 = 0 

with 

X s - p 0 cos * 0 
y s - p 0 sin * 0 

6^ = a b - p q (b cosi> 0 cos v + a sin 4 sin v) 

Bo = b sina cos v - a cosa sin v 

c. S S 

& 3 = b cos<i> s cos v + a sin<f> s sin v 
2 2 

e 4 = (b -a )cos v sin v + p Q (a cos<f> 0 sin v - b sin<|> 0 cos v) 

Now using the substitution 
a = e jv 

or in other words 
a 2 +l 

cos v = .-25- 
and 

sin v ^ 
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then 


1 2 
= 2a {2abct " p o^ b C0S V j a sir V“ + ( b c0S< J> 0 + j a sin«j> 0 )]} 

1 2 

6 2 = 2a ^ b sin ^ s +J ' a cos 4> s )“ +(b sin<j> s - ja coscj> s )] 

1 2 

B 3 = 23T ^ b cos< J’ s "J a sin< J> s )a +(b cos<j> s +ja sin<f> $ )] 

64 = [J'(a 2 -b 2 )a 4 -2p 0 (b sin^+ja cos^Ja 3 

4a 

2 2 

- 2p o (b sin<j> 0 -ja cos4> 0 )a - j(a -b )] . 

Multiplying the equations together, collecting terms and setting 
the equation to zero results in a sixth order polynomial in a 
given by 

Cga^+Cga b +C^a^+C2a 3+ C2a 3+ C-|a^+C o = 0 

where 

2 2 

c 6 = (a -b )(a si n4> s +j b cos<j> s ) 

C 5 = " 2 P 0 [( a2+b2 ) s -' n (V^s )+J ' ab C0S (V^s )] 

C4 = a(5b 2 -a 2 )sin<}> s + j b(5a 2 -b 2 )cos<j> s 
C 3 = 4p (a 2 -b 2 ) sin( 4> 0 +4> s ) 
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C 2 = (complex conjugate of C^) 




This can be solved using any standard numerical method for the 
solution of polynomials with complex coefficients. Of course, 
there will be six roots to this equation, but only one will be 
physically realizable. This one is found by plugging the v's 
from 

v . tan -l (iSMinar^ 


into the equation for I and finding the one which gives the 
minimum distance |I|. The parameter z e can now be found from 
Equation (66b), that is 


z s 


( Vx * n vV d z 
<"x d x + n yV 


( 67 ) 


The method for finding the reflection point discussed above can 
be considered an exact solution. There will be some errors involved 
in practice due to the numerical solution of the polynomial, but these 
will not be accumulative since the polynomial must be resolved for 
every scatter direction. This method is time consuming, however, 
because it must be resolved every time a new pattern direction is 
considered. Also, more complicated interaction problems involving 
reflections, as discussed in the next chapter, will not have such 
a simple solution. It is for these reasons that an approximate 
method for solving for the reflected field is desirable. 
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In most radiation pattern calculations, the scatter direction 
is normally incremented by only a small amount before another 
field point is desired. This suggests the possibility of using an 
incremental scheme to find the reflection point on a curved surface. 
This method is developed based on the first term in a Taylor expansion 
of the reflection equation 


f(v "»s ) = <n x I x +n yV (n xV"y d x ) + < n x d x +n y d y )(n xWx ) = 


( 68 ) 


that is 


Vi ■ f J + df : 


where 


9f • 3f j 

df j = s/ 9V + ^ 


Since 


f w 


Vj + SV 

V^s 

this implies that 


= 0 


or 


f. + df, 
J J 


V j 


= 0 


3f • 

f . + — -J- 

J 9V 


6V + 


Hi 

94> c 


5<j). 


= 0 


so that 
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Sv = 


-f . 
J 


■57— <$<}> 
3 »s s 


3 V 


This says that if the reflection point for an observation angle of 
4> s is at the elliptic angle Vj on the cylinder surface, then the 
reflection point for an angle of <J> S +I $<I> S is given by 

Vi c v j +sv - 

From Equation (68), it is seen that 


3f. / 3d 

- 1 ■ (M +n„IJ n. 


3d 


3<f> 


XX y y' l x 34 


- n 


y 3 <i> 


f 3d x 3d \ 

+ ("xC + n y ) (n x‘y ■ n y’x> 


and 


3f . / 3n. 


— = — 2. I + n 
3V V 3V X "x 


+ (n d +n d ) 

' v v \i \t • 


3l 

3n 

X 

+ I 

3V 

3 V 

3n 

X 

I + n 

3V 

y x 


31. 


!Lt . , ik 

3V 3 V X y 3V 


/ 9n y 3n 

+ (n I +n I )( d - — d 

1 x x yy'Vsv y 3v x„ 


3n. 


3n 


V^3v d x + 3v d yy ^ n x*y _n y*x^ 
and from Equations (54, 63) and (65a) 
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3n. 


av 


x , 

=-b sin v 


jr = a cos v 


'5 

av 


= -a sin v 


ai 


av 


— = b cos v 


3d 

= -sin<J> 
a<p- s 


3^ = C0S V 


Now Vj + i may be solved for, and an error computed from 


e = n . 1 1 + n • d 

j+1 


V’ 

* s +« s 


If the error is not within a given bound, the number of interval 
steps 64> s may be halved and the total number of steps doubled so 
that the new scatter direction is the same but the increments have 
doubled. This process may be continued until the error bound is 
reached or the roundoff error in the computer has been reached. 

The z component of the field may be found from Equation (67) 
as before. 
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The assumption in the above method is that the previous 
reflection point is known accurately. Of course this is not true 
in general, since the previous reflection point is computed 
using the same method from its previous point. This means that 
the process can accumulate error, but this is compensated for 
by allowing the step size to vary and the fact that the error 
is computed against the laws of reflection that must be obeyed. 

In other words, the errors accumulated in calculating the previous 
point (fj) are included in the calculation of the next point (fj+i). 
Thus, as errors grow the interval decreases and tends to correct 
the total result at (fj + -|). 

A starting location for the first computation needs to be 
provided for the incremental method. This is most easily ob- 
tained from the fact that the reflection point is known at the 
shadow boundary, since this is simply the point of tangency from 
the source to the cylinder. The two points of tangency v-ji and 
vi2 were previously found in Equation (60). The scatter direction 
at the start will be the tangent vectors -| and t^ 2 » respectively. 
The one closer to the desired direction cT can be chosen to start the 
incremental method. 


The reflected field from the elliptic cylinder is found using 
Equation (2). The geometry is illustrated in Figure 40. In this 
case the observation point is in the far field, Using the 

usual far-field approximation and referring the phase to the origin 


f(e. 




r 


(Q r ) 




r r 
jP£ 


J'kR -d 


- j kd 


(69) 


where R r is the reflection point and R is the reflection coefficient 
given by Equation (3). The principal radii of the reflected fields 
for spherical wave incidence are given in Equation (10), where for 
an elliptic cylinder R 2 -**>. This gives 

p£ = s' (70a) 


s' R, cose 1 

L -j o— • (70b) 

R-|Cose +2s'sin e 2 
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where 


(71a) 


(b 2 cos 2 v +a 2 sin 2 v ) 3/2 

R 1 = Hi 1 


(71b) 


and 


cosQ = n • d 


2 2 2 2 1 
sin e 2 = sin w + cos to cos 6 


(71c) 


(71 d) 


and 


cosu)y, — -I • e 2 


sinu> r = -I • e-| 


e, - e 2 X " = e lx x + e ly y 


A A 

e 2 = z 


(71e) 
( 71 f ) 
(71 g) 
(71 h) 



Figure 40--Reflected field geometry for elliptic cylinder. 
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The vectors that are used in the reflection coefficient are given 

by 


e x = sin(o) r -90°)e lx x + sin(o) r -90°)e ly y + cos(uy90°)z, (71 i ) 

A • A A 

e il = I X e j. ’ (71 j) 

and 

A /V A 

ej = d x e x . (71k) 

The incident field ^(Qd) is given as in Equation (1), with the in- 
cident vector I defining the incident angles. 

Reflections can also occur off the flat surface of the end caps. 
The normals to the end caps are defined as 

~+ + * 

ri = -cose x + sin9„ z (72a) 

ec c c v 


and 


_ A 

n„„ = cose x - sine z, 
ec c c 


(72b) 


where the positive sign refers to the most positive end cap with its 
origin at z = zj and the negative sign refers to the most negative 
end cap with its origin at z=z“. It is possible for reflections to 
occur if 


where 



If the above inequal ity_is true, then an image can be found and 
an intersection point (R^) on the plane of the end cap can be 
found from Equations (31 and 33) that were used for the flat 
plates in Chapter III. The test to see if the point of inter- 
section lies in the confines of the end cap, however, is slightly 
different. First, the distance from the origin of the end cap to 
the intersection point of the end cap in the plane of the end cap 
is defined, that is 
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If 



P t < min(a csce c , b) 

then a reflection will occur, or if 
P tc > max(a csc0 c , b) 
then a reflection will not occur. If 

min(a csc9 c ,b) < p t < max(a csc0 c , b) 
then another test is constructed as follows. Let 

R t ■ R tx i + V + R tz* 

then 

-,/a csc 0 R t • y 
v = tan - ( 5 

V b R t • x c 



where 


*c -< 


y x n 


ec 


vec 


The radius of the end cap at the angle v c is then given by 

- 2 5 2 

p = / a esc e cos v + b sin v_ . 
c v c c c 

Thus if 

p tc < p c 

reflections occur. 
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The reflected field for both end caps, if they exist, may be 
found by the same procedure used for the flat plates in Chapter 
III, Equations (35-39). The total reflected field from the 
end caps can be written as 


V' . (e ,d> ) = E r+ . (e ,d> ) + E r " . (e ,<j> ), (73) 

end cap' s Y s ; end cap' s w s' end cap s v s' ' ’ 


where the plus sign refers to the most positive end cap and the 
negative sign to the most negative end cap. Usually both end caps 
will not contribute to the field at the same time. 

Curved surface diffraction is used for aspect angles in the deep 
shadow region. The geometry is illustrated in Figure 41. This 
diffraction coefficient is discussed in Section G of Chapter II. 


SOURCE 



Figure 41 --Curved surface diffracted field geometry. 
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When the scatter direction is in the far-field, Equation (23) is 
given by 


eV,.,*,.) = E i (Q 1 ) 


.. jk R.-d 

[b-,b 2 u s + n 1 n 2 u h ]Jp e 


e -Jkd 

d 


(74) 


where R<j is the point at which the ray detaches from the surface 
(Q 2 ) and sheds in the scatter direction. The functions u s and 
u^ are defined in Equation (24) and p is given in Equation (25b) 
for a cylindrical surface. 

The point Q-| is the incident point and it is found from 

Equation (60). There are two incident points for a closed convex 

surface such that energy propagates in a clockwise (cw) and 

counterclockwise (ccw) sense around the surface so 
r 


il ’ 

ccw 

i2’ 

cw . 


The surface ray follows a geodesic path along the surface, which 
is a helical path for a cylinder. The point at which the energy 
sheds off the cylinder is given by 


tan 


,-| (- b cos$ s ' 


a sin<j> 


v f 


tan 


/ b cos<f> s 




-a si n 4 > s ) ’ 


ccw 


cw 


If one defines the geodesic starting direction by the angle a = 
180°-e s , then 


z 


- C0Sa s 

|sina s i 


f v f r~2 2 2 2 

1 vj a sin v + b cos v dv 



(geodesic equation) . 


110 



This equation can be quickly evaluated using numerical techniques. 
The important parameters of this problem are listed below: 


t = 


sina 


f f 1-2- 

Ja sin 


2 2 2 

v+b cos v dv (arc length) 


v- 


; _ -a sin v x + b cos v y 

1 r-? p 2 2 

yj a sin v + b cos v 


e 2 2 


f (curval inear coordinates) 



sin Vl 


A 

COSa^e-, 
s 2 


v=v 


i 


v=v 


f 


(unit tangent vector) 


1,2 


_ b cos vx + a sin vy 


/ 


2 2 2 2 
a sin v + b cos v 


v=v. 


v = v , 


(unit normal vector) 


b l ,2 = h ,2 x n l ,2 


cosa s e-| + sina s e 2 


(unit binomial 


v=v.j vector) 


v = v , 
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(longitudinal radius of curvature) 


(a 2 sin 2 v + b 2 cos 2 v) 3/2 

, . 2 
a b sin a 

s 


p t 


(a 2 sin 2 v + b 2 cos 2 v) 3|/2 
ab | sin 2 (a g - 90°) | 


(transverse radius of curvature) 


R , = a cos v f x+ b sin v r y + (z +|l I cote + t cose )z . 
d f f s 1 xy 1 s s 


Using the above relations. Equations (74) can be employed to find the 
surface diffracted fields from the clockwise (E° ) and counter- 
clockwise (Eccw) propagating rays. The total surface diffracted 
field from the cylinder, away from the shadow boundaries, is 
given by 


= f c„<V*s> + Ov*s> 


(75) 


As stated previously, both the reflected field and the surface 
diffracted fields for the elliptic cylinder are inaccurate close to 
the shadow boundaries. The transition field solution of Section H 
in Chapter II for a curved surface is needed in this region. The 
transition field for a cylindrical structure with the scatter 
direction in the far field is given in Equation (26) as 


E t * E^ (Q) • [u 1 u 2 G s + Vl v 2 G h ] 


(76) 


where r 



for the lit region 
for the shadow region. 


In the lit region, the functions (G and G, ) are given by 
Equation (27), as 
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and for the shadow region Equation (28) gives 
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(77b) 


The necessary equations to define these functions have been given in 
Chapter II, Section H and this chapter, where the parameters in the 
reflected field section are used for the lit region and the param- 
eters in the surface diffracted field section are used in the 
shadow region. 

The transition field solution is accurate at and around the 
shadow boundaries. As the observation point moves away from the 
shadow boundary either into the lit region or deep shadow region, 
the geometrical optics fields and surface diffracted fields, 
respectively, present a more accurate solution. At some point 
the solutions are very close together and either solution is valid. 
This is illustrated in Figure 42 for an elliptically capped thick 
wall (a=1.59X). The dipole is on the x-axis, 4.76 from the 
origin, and parallel to the z-axis. This indicates that the 
transition region is rather arbitrarily defined in terms of the 
accuracy of the solutions. A progressive blending function can 
be constructed, so that the total solution smoothly goes from 
one type of solution to the next over a small region in the 
area of common validity. 
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WALL 


Figure 42a--Cyl indrically capped wall configuration excited by 
a dipole source parallel to z-axis. 


For a closed convex surface, there are two shadow boundaries 
as discussed before so there will be a transition field associated 
with the top or counterclockwise term (E£ c ) and the bottom or 
clockwise term (e£ w ). The total transition field, where it exists, 
is then given by 


,<f> ) = E* (9 ,(f> ) + E* (9 >( f> ). 

s s' cw s Y s ccw' s *s 


For an infinitely long cylinder, where the end caps need not 
be considered, the total field is given by 


E(9 S ,<)> S ) = E 1 +f r +E t + E d 
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TRANSITION 



figure 42b. Eg radiation pattern of electric dipole source 
surrounding cylindrically capped wall showing 
the effects of the three different solutions 
in their respective regions. 
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TRANSITION 



Figure 42c--E ( j ) radiation pattern of magnetic dipole source 
surrounding cylindrical ly capped wall showing 
the effects of the three different solutions 
in their respective regions. 
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using Equations (29,69,75,78). The accuracy of this result is 
confirmed by numerical comparison with the results obtained via 
an exact eigenfunction solution for a circular cylinder [64] in 
Figures 43 for the acoustic soft and hard boundary conditions, 
respectively. The results for the GTD solution are also com- 
pared against a moment method (MM) solution for the elliptic 
cylinder case in Figure 44 for the hard boundary condition and 
two different source positions. The result in Figure 45 compares 
the GTD solution with the MM solution for a conical scan angle of 
0 = 45°. In all the above comparisons the results show excellent 
agreement. 



Figure 43a--Circular cylinder configuration excited by a 
dipole source parallel to z-axis. 
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SOFT b.C. 


Figure 43b--Ee radiation pattern of electric dipole surrounding 
circular cylinder. Here, GTD implies the GTD + 
transmitter solution. 
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GTD 

90* « ® • • EIGENFUNCTION 



HARD b.C. 


Figure 43c— E a radiation pattern of magnetic dipole surrounding 
circular cylinder. Here GTD implies the GTD + 
transition solution. 
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SOURCE 



a * 2\ b * I X , />' * 4 X 


Figure 44a--Elliptic cylinder configuration excited by a 
magnetic source parallel to the z-axis. 
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/>'«= 4 X 
HARD b.C. 


Figure 44b--E<j> radiation pattern surrounding the elliptic 
cylinder. Here GTD implies GTD + transition 
solution. 
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90 ° 



Figure 44c — Etj, radiation pattern surrounding the elliptic 
cylinder. Here GTD implies GTD + transition 
solution. 


122 






ELLIPTIC CYLINDER 



Figure 45a--3-D scattering configuration involving an infinitely 
long elliptic cylinder. 
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OUR SOLUTION 


90° XXX MOMENT METHOD 



Figure 45b--Conical scan radiation pattern for the 3-D 

scattering configuration. The pattern of the 
<j> component of the electric field is shown. 
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The flat end cap and elliptic cylinder junction can be 
represented by a curved wedge such that one can apply the 
Geometrical Theory of Diffraction as discussed in Chapter II, 
Section E. The end caps are defined by their normals as given 
in Equations (72a and b) , and their origin position on the z- 
axis. The edge of the wedges are given by 

- a + . , ~ . . + 
e = - a sin v x + b cos v y - a sin v cot e c z 


for the most positive end cap and 

e^ = a sin v x - b cos v y + a sin v cot e~ z 

J c 


for the most negative end cap. The diffraction by a curved edge 
obeys the same laws of diffraction as the straight edge in that 
the half-angle for a cone of diffracted rays is equal to the 
angle of incidence. This can be expressed as 

A A A A 

e(v) • I(v) = e(v) • d . 


The location of the diffraction points is more complicated, how- 
ever, since the edge vector is a function of the parameter v. It 
can be shown that for a circular or elliptic edge outside the 
caustic regions there are in general four points of diffraction on 
the edge of a near field source and a given scatter direction. 

The diffraction points can be found in many ways. An exact 
solution similar to that obtained for the reflection off an 
elliptic cylinder can be found which in this case is an eighth 
order polynomial. An incremental solution could also be found 
for this problem, but this is complicated by the fact that three 
of the diffraction points can merge into one in certain scatter 
directions. A search method is another possibility and is used 
here. The geometry is illustrated in Figure 46. 

The search method is based on minimizing the error given by 


e(v) = e{v) • I(v) - e(v) • d 

where 

T(v) = (a cos v - x $ )x + (b sin v - y $ )y + (a cos v cot e c 

+ z r - zjz 
c s 

as v is incremented from 0° to 360°. A minimum is found using 
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SOURCE 

LOCATION 



DIFFRACTION POINTS 
ON RIM 


E *2 


Figure 46- -Curved wedge diffraction points on rim of end cap 
of finite elliptic cylinder. 


= e( v j_i) + e( v j) + e(Vj +1 )» j = 0 to 361. 
If 


Ujl i Uj.,1 and |Cj| < |c j+ ,|. 

then vj is the closest point to a minimum. An interpolation 
scheme J is used to find the actual minimum based on 

V d = V. + AV 

where 


AV = s' (e iJL.-J ^ I_ __ 

e*Al+Ae.I-(I-Al)(e*<?)-s 1 (Ae’d) 

and 

s' = |T| 

A A A 

aI = -a sin v x + b cos v y - a cos v cote c z 


- 1 A A A 

Ae” = ±(-a cos v x - b sin v y - a cos v cote” z). 


This process is continued until four points are found or the whole 
edge has been searched. 


Now that the points of diffraction are known the diffracted 
field can be determined. In the far field, the diffracted field 
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The diffraction coefficient for the curved wedge are given in 
Equation (19). The distance parameter associated with the in- 
cident field with s-*» is given by. Equation (20a), 

L i _ <^2 sl ' n2 B 0 
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For spherical wave incidence 



so that L' is given by 
L 1 = s' sin^B Q . 


The distance parameter associated with the reflected field with s-**> 
is given by. Equation (20b), 

r r -2 
,r = P 1 p 2 S1n e o 
r 


v* 

where p is given in Equation (21). For the 0 surface, which is the 
flat ena cap, the parameters are simply given by 


p 


ro _ 
1 




= s' 



sin^B 


o * 


For the n surface, which is the curved surface, the parameters 
pC n and p£ n are defined in Equations (70a) and (70b), respectively. 
The parameter p£ n is given by 



2(n-n e )(n-i) 
a e si " 2e o 


where n is the unit normal to the cylinder at the diffraction 
point, the normal to the edge is 
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n = b cos Vasins 
e d c 


x + 


a sin 

iTnT 

c 


y + b cos v^ cos 0 c z, 


= (a 2 sin 2 v+b 2 sin 2 e c cos 2 v) 3 ^ 2 

o 

ab sin 0 C 





= a cos Vj x + b sin y + (a cos v^ cot e c + z c )z. 


The caustic distance p is given by. Equation (18), 


1 

p 


1 _ 

s' 


A A A 

V (I - d) 

r 2 

a„ sin 6 
e o 


As for the straight edge a ray fixed coordinate system is used 
centered about the diffraction point as shown in Figure 47. The 
necessary equations^are given in Chapter III, where here e cor- 
responds to e m and n ec to n|pi a t e . The wedge angle number for 
the curved edge is given by 


n w = 1 + 


-1 ~ 
cos (n 


n ec> 


v=v. 


If <j>* > n^ (180°) or ip > n w (180°), then the incident or dif- 
fracted field, respectively, is shadowed and does not need to be 
found. In other words, it is assumed that only diffraction 
points that have an incident field direct from the source is to 
be found. It is possible to have a diffraction from the curved 
edge that results from a field that is incident on the edge after 
it has propagated along the cylinder as a surface diffracted field. 
This type of field is neglected in this study. The diffracted field 
from the jth diffraction point on the ith end cap is then given by 
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E?j(V*s> = E ^o + E l+- 


As is mentioned above, it is possible for three of the dif- 
fraction points to merge into one for certain aspect angles. This 
represents a pseudo-caustic effect resulting from the disappearance 
of isolated ray contributions. The GTD fails to describe the 
transition from the four isolated diffracted ray result to the 
pseudo-caustic result and hence needs modification. This is 
presently being studied [65] through the use of equivalent rim 
currents which indirectly employ the GTD edge diffraction coeffi- 
cients. The equivalent edge currents are then incorporated in the 
radiation integrals which yield a smooth behavior for the field 
across these pseudo-caustic regions. At the present time, this 
has not been included here. 



Figure 47 — Edge coordinate system at curved wedge diffraction 
point. 
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The total diffracted field from the end cap rims, if they 
exist, can be written as 


r 1 


2 4 


end cap (e s *^s } = ^ ^ij * 


(79) 


The total scattered field from a finite elliptic cylinder is 
now given by the sum of all of its individual contributors in 
their appropriated regions as 

■ f1 + ^ + ^nd cap + + ft + ^end cap. < 3 °> 


The description of the individual scattering contributors from 
the cylinder are summarized in Table III. 


TABLE III 

SCATTERED FIELDS INCLUDED' IN THE FINITE CYLINDER SOLUTION 



Symbol 

Description 

Figures 

Equation 

1) 

r 


Reflected field off the 
cyl inder 

40 

69 

2) 

r 1 


Diffracted field off the 
cylinder 

41 

75 

3) 



Transition field off the 
cyl inder 

40 & 41 

78 

4) 

f r 

end 

cap 

Reflected field off the 
end caps 

- 

73 

5) 

E end 

cap 

Diffracted field off the 
edges of the end caps 

47 

79 
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To illustrate the validity of this near field solution, it is 
compared with measured results for a finite circular cylinder given 
in Reference [23]. The geometry is illustrated in Figure 48. For 
a horizontal dipole parallel to the cylinder axis, the present 
result is compared with a measured result as shown in Figure 49. 
Note that the two results are almost identical except about the 
caustic region where the analytic solution does not accurately 
describe the field. The same is true for a vertical dipole 
(parallel to the y-axis) result shown in Figure 50. The caustic 
region is not extensive for most cases. Leaving this caustic in 
the solution, does not strongly effect the results in practice 
as will be seen in the next chapter in which the aircraft problem 
is studied. 



Figure 48--Short dipole illuminating a finite cylinder. 
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Figure 49--Radiation pattern (Eg) of a horizontal dipole (parallel 
to z-axis) in the x-z plane. 
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Figure 50--Radiation pattern (E^) of a vertical dipole (parallel 
to y-axis) in the x-z plane. 
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CHAPTER V 

ANALYSIS OF WING MOUNTED ANTENNA RADIATION PATTERNS 
A. Introduction 


The previous chapters have illustrated how the Geometrical 
Theory of Diffraction is used to study the scattering by basic 
structures (flat plates and finite cylinders). In this chapter 
these solutions are extended to determine the radiation patterns 
of antennas mounted on the wings of aircraft. The analytic 
aircraft model consists of a finite elliptic cylinder representing 
the fuselage with flat or bent plate wings and horizontal and 
vertical stabilizers. Small plates and elliptic cylinders can be 
used to represent engines and or external stores if necessary. 

The basic aircraft shape to be modelled is illustrated in Figure 
51. This model is versatile enough to simulate the major scattering 
contributors of a wide variety of aircraft. The solution is also 
formulated in such a way that an arbitrary antenna can be considered 
simply by integrating the equivalent aperture currents. The lower 
frequency limit of the solution is dictated by the ray optics format 
which requires that the various scattering centers be no closer than 
approximately a half wavelength with the overall aircraft being large 
in terms of the wavelength. The upper frequency limit is dictated 
by the model representation of the actual aircraft considered. 

B. Wing Mounted Antenna Analysis 

In general, an antenna mounted near a number of flat plates and 
an elliptic cylinder is even a more complicated scattering problem 
than that of a number of flat plates alone as discussed in Chapter 
III. Again, it is necessary to make a few good engineering assump- 
tions about the problem in order to make it tractable. First it is 
assumed, at present, that the antenna is mounted on or close to the 
wing structure, and it is at least a small distance (greater than a 
quarter wavelength) from the fuselage. This is done merely to reduce 
the number of scattering mechanisms to be considered. An antenna 
mounted on the fuselage needs a different treatment as was originally 
shown in Reference [23]. When the antenna is mounted very close to 
the fuselage neither solution is applicable in its present form. 

The antenna is also assumed, at present, to be mounted away from 
engines and external stores. If the antenna is mounted sufficiently 
far from these objects , (i .e. , the antenna is mounted On the top 
portion of the wings) their effects will be negligible. If, however, 
the antenna is mounted on the lower portion of the wing and is close to 
an obstacle, the scattering effects of that object can be approximated 
by using small flat plates or elliptic cylinders. This is not done in 
the present analysis. It should also be possible to easily extend 
the present solution to antennas mounted on the vertical or horizontal 
stabilizers. 


135 






Figure 51--IT]ustration of general aircraft model for wing 
mounted antennas. 


136 


The scattering mechanisms considered for this model include only 
first and second order interaction terms with second order being mecha- 
nisms like reflected-diffracted or diffracted- reflected fields. All 
the scattering terms considered for the scattering from a number of 
plates in Chapter III, Table II, and for the finite elliptic cylinder 
in Chapter IV, Table III, have been included in this aircraft model. 

They will not be repeated again. The additional interaction terms 
that have been found to be important for the aircraft model are shown 
in Table IV. Subscripts have been added in the notation of the field 
to distinguish between interactions between the plate (pi.) and the 
cylinder (cyl.). Note it is assumed that the wing on the opposite 
side of the structure from the antenna does not contribute to the 
scattered field. This is a valid assumption in that the surface ray 
must propagate sufficiently far around the fuselage in order to 
illuminate the other wing which implies the surface ray attenuates 
a great deal. Thus, the field illuminating the opposite wing is 
very weak as well as its scattered field. Also, the wing is as- 
sumed to be mounted in the y-z plane, as illustrated in Figure 51, 
that is at the symmetry point of the cylinder so that image theory 
can be used to model the wing-cyl inder junction. This can be modi- 
fied by the use of curved wedge diffraction, if the wing must be 
mounted at an angle. However, this does not appear to be necessary 
for a practical representation of real aircraft configurations. 

As mentioned before, it is important to shadow the fields properly. 
This can be done by checking the incident ray paths and scattered 
ray paths for all the interaction mechanisms using both the plate 
and cylinder shadowing procedures outlined previously. Only those 
fields that exist should be included in the total field, the shadowed 
fields must be set to zero. This will help insure that significant 
terms that are large enough to produce visible discontinuities in the 
radiation pattern when neglected will be detectable. These higher 
order terms can then be added to the pattern if they are practical 
to obtain and are of engineering value in the final result. 

Each of the scattered field terms that are to be included for 
the aircraft model, as listed above, involve a scattering mechanism 
in which one part of the structure acts as the source of the incident 
field that illuminates a different part of the structure. The first 
three scattering mechanisms listed above are simple to determine 
because they involve reflections off the flat plates such that 
image theory can be used. As a consequence, the various mechanisms 
can be viewed as point sources radiating spherical waves originating 
from their images. This means that the analyses carried out in the 
previous chapters for the individual scattering terms are applicable 
here. The procedure for obtaining these fields need not be re- 
peated here. The last two scattered field terms in the list, 
however, are -compl icated by the fact that after the first interaction 
with the structure, the incident field for the second interaction 
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can not be viewed as a point source with a spherical wavefront. 

These terms must be treated as a source that radiates a wave with two 
radii of curvature. Consequently the new problems associated with 
these two scattering mechanisms are considered in great detail. 


A 



Figure 52--I1 1 ustration of reflected field off the cylinder after 
reflection by a plate. 


The first scattering mechanism to be considered is the reflected 
field of the cylinder after reflection from the jth plate. This 
is illustrated in Figure 52. The reflected field of the cylinder 
is given by Equation (69), where the origin of the incident field 
is the image of the source (R-j) in the jth plate. The incident 
field then is the reflected field of the jth plate in the direction 
from the image to the reflectiQn point of the cylinder and is 
given by Equation (34), with ^(Qe) = E^e-j ,<j>.j ) . The application 
of Equations (37-39) to find the reflected fields of the jth plate 
and the applications of Equations (70) and (71) to find the,, 
reflected fields off the cylinder in the scatter direction d is a 
straightforward extension of the above equations and need not be 
discussed further here. The total field associated with this 
scattering term is given by 


-JT 

Epl-cyl ^ 9 s’^s^ 


E j,cyl ^s’^ ’ 


(81) 


__y» 

where E . , (e ,<j> ) = 0 if it is shadowed in any way. 

J 9 cy i s s 
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Figure 53--I1 lustration of transition field off the cylinder after 
reflection by a plate. 


The scattering mechanism to be considered next is the transition 
field scattered by the cylinder after a reflection by the jth plate. 
This is illustrated in Figure 53. The transition field of the 
cylinder is given by Equation (76), and again the image of the 
source in the jth plate can be viewed at the origin of the incident 
field. The incident field for this term is then given by Equation 
(34) which is the reflected field from the jth plate in the direction 
from the image in the plate to the cylinder scatter point. The 
cylinder scatter point in the lit region is the point of reflection, 
and in the shadow region it is the tangent point. The necessary 
equations for the reflected field from the jth plate are given in 
Equations (37-39) and the necessary equations for the transition 
field given in Equations (77a and b). Note that the shadowing 
procedures determine if the transition field due to both the top 
part of the cylinder and the bottom part of the cylinder need to 
be included. The total field associated with this scattering 
term is given by 


Ft 


p p 

•prt A v -?rt 


, (e ) = y F. L +J . 
pl-cyl' s ,M V j,cw j ,ccw 


(82) 
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The scatter terms associated with a reflection off the jth 
plate after being reflected by the cylinder are illustrated in 
Figure 54. The reflected field of the cylinder is found in the 
direction obtained by imaging the scatter direction into the plate 
using, 

d r = d - 2 (n-d)n . 


A 



Figure 54-- 1 1 lustration of reflected field off, a plate after re- 
flection by the cylinder. 


This reflected field, given by Equations (69), is the incident 
field on the plate such that 


The reflected field off the plate is then computed using Equations 
(37-39) in the scatter direction. The field appears to be coming 
from the images in the plate of the two caustics associated with 
the reflected field off the cylinder. The total field for this 
scattering term is given by 


-rr 


'cyl-pr s ,y s 


(0 e »O = I 


j = l 


^cyl ,j ^ 9 s’*s^ 


(83) 


The scattering mechanism that corresponds to a diffracted field 
from an edge after a reflection off the cylinder is illustrated in 
Figure 55. In order to find this field in a given scatter direction 
d, the ray path that the field travels must be determined. This 
path is governed by the laws of reflection off the cylinder surface 
and the laws of diffraction along the edge. 

The law of reflection on the cylinder can be written as given 
in Equation (64) such that 


(nxl)(n-r) + (n-I)(nxr) = 0. (84) 


In this case the ray designated r is not known a priori. It is de- 
pendent on the diffraction point on the edge. The diffraction 
point, in turn, is also dependent on the reflection point. This 
means that the above equation is non-linear. At present, there 
does not appear to be a way to solve this problem in an exact 
manner. The incremental method, based on a first order Taylor 
series approximation, as developed for the reflection point off 
an elliptic cylinder in Chapter IV may be employed here. 

For a cylinder the laws of reflection as governed by Equation 
(84), reduce to the two equations 


f(v,z,e s ,* s ) = (n x I/n y I y )(n x r y -Y x ) + (n x r x+ n y r y )(n x I y -n y I x )=0. 


and 


9 (v,z,e s ,* s ) = (n x I x +n y I y )V(n x r x+ n y r y )I z =0. 


(85a) 

(85b) 
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TOP VIEW 


Figure 55-- 1 1 lustration of diffracted field off an edge after re- 
flection by the cylinder. 
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The normal to the cylinder is given by 
n(v) = n x (v)x + n y (v)y 
or 

n(v) = b cos v x + a sin v y . 

The incident vector on the cylinder is given by 
t(v,z) = I x (v,z)x + Iy(v ,z)y + I z ( v ,z)z , 

t(v,z) = R r (v,z) - R s , 
or 


f(v,z) = (a cos v - x s )x + (b sin v-y $ )y + (z-z $ 

where R r is the reflection point and R s is the source 
The reflected ray is given by the vector 

r(v,z,9 s ,<f> s ) = r x x + r^y + r z z 
or 

F(v,z,e s ,* s ) = R d - R r 

where Rd is the diffraction point on the mth edge and 

Rj = C m + t Av,z)e m + s (v,z) cots (e,. , 4 >^. )e m 
d m p v ' m p ' o s r s m 

and 

t_ = (R -C_) • e m , 
p ' r m m 

s = IR”, - t„ e m - LI , 
p 1 r pm m 1 


z 

location. 


is given by 
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cots o 



Since r is dependent on v and z. Equations (85a and b) are not 
separable and must be solved together. Expanding these two 
equations in terms of a first order Taylor series gives 


f j + l ’ f j + w 1 5V + ^ 62 + 5e s + ^ + 0(52 > 


ag. ag. ag. ag. ? 

g j + l - g j + -iw 4V + Si 1 SZ + 3^ 6e s + + °< 6 > 


Since from the laws of reflection 


Vi - 0 

g j + i = 0 



This two by two matrix is easily solved to give 


Sv 




9f. ag. 


az ae. 




+ 




(86a) 
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(86b) 
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The partial derivatives are given by 
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with 
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1 x 3z y dz J z 'xx y y 3z 
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^ n x I x +n y I y^ 33 


3r 


n 


x 33 


3r 

+ n -X- 
y 3B 


with 3 = 0 s , 4> where 



-b sin v 


i = 

3v 


COS V 


i 1 * 

3v 


-a sin v 


31 

= 

3v 


b cos v 



= 1 



3a 3a 3a 


with a = v, z 



3v 


-a 


sin v x + b cos v y 
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5 

3z 


3R d at 


9s p 

aa a« e m + a a cot e 0 e m 


I!e. = . e 

3a 3a m 


3S. 


|3 _ 3a 

3a 


3R r /n 7 T x ( sR r e )( R - C ) • e 

— • m A r nr m 


<V C m> H 3a 


/ 


(R -C ) • (R — C ) - [ (R -C ) • e V 
r m r m r m m 


a? _ 9 (c°te 0 ) , 

33 s p 33 e m 


with 3 = e s , <J j s and 


3(cot3 Q ) 

33 


3d 

■ 33 ’ e m 

[Md-e /] 3 ' 2 


These equations can then be used to find 6v and 6z changes 
associated with the reflection point on the cylinder surface due to 
60s and 6$ s changes in the scatter direction. The new reflection 
point is found from the previous reflection point using 



+ 6v 


and 


z 


j+1 


= Z. + 6Z . 
J 
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Once the reflection point on the cylinder is known the diffraction 
point on the edge can, also, be determined. 

Errors can be computed to see how accurately the values are de 
termined using 


e 


C 


n 


v 1 

2 j +1 


+ n 


r 


Vr 

0 s+ 66 


Z j+1 

s’*s +6 +s 


and 



v j+l ,z j+l 

e s +60 s**s +6 *s 



d 


e +69 
s 

4> s +5<f> 


s 

s 


If the errors are not within a predetermined bound the interval 
steps 60 s and 6<j> s can be halved and the procedure repeated until 
the desired scatter direction is obtained and the desired error 
bounds are achieved. 

A known starting location needs to be obtained for the first 
computation of the incremental method. This can be found in this 
instance by noting that the corner points of the edges are known 
and in fact they provide definite bounds on the permissible range 
for the reflection-diffraction path. With the corner points known, 
as shown in Figure 56, Equations (85a and b) can be solved by means 
of the polynomial solution technique derived in Chapter IV. 


If the mth corner is given by 


C = 
m 


V + y c y + z c z 


then 


r = (x c -a cos v)x + (y c ~b sin v)y + (z c ~z)z 


with n and I defined above. These values are substituted into 
Equation (85a) along with the fact that 
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Figure 56--Illustration of starting point for finding diffraction 
point on an edge and the reflection point on the 
cylinder for the reflection-diffraction problem. 


150 


The equation is multiplied thru and like terms in a are collected so 
that a sixth order polynomial is obtained given by 

6 5 4 3 2 1 

CgCt + Cgd + C^a + C^a + Z^l + C-jCi + C Q =? 0 

where 

C g = (a 2 -b 2 )[a(y c +y s )+jb(x c +x s )] 

C 5 = -2( a2 +b 2 )(x s y c +ys x c)-j 4ab[(a 2 -b 2 )+(x s x c -y s y c )] 

C 4 = a(5b 2 -a 2 )(y c +y s ) + jb(5a 2 -b 2 ) (x c +x s ) 

C 3 = 4(a 2 -b 2 )(x s y c+ y s x c )- 

C 2 = (complex conjugate of C^) 

C 1 = C 5 

C 0 * C 6 ' 

The six roots of the polynomial correspond to v in that 



The physically realizable value which is the true reflection point 
is found by finding the v that minimizes the distance given by 
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+ M » 


which is necessary to satisfy Fermat's principle. Using this 
parameter (v r ) representing the reflection point in the x-y plane, 
the z-coordinate can be found from Equation (85b) to be 


2 KyyVvKyVy) 2 .: 

r (n x r x + y y >-<Vx +n yV 


With the reflection point on the surface of the cylinder known 
for a given corner on the mth edge, the vector r represents a bound 
on the possible incident rays on the edge. The angle of the dif- 
fracted cone of rays from the mth edge at the corner is given by 


One of the rays defined by this cone angle can be used as the 
starting scatter direction for the incremental method. It is 
defined by 



(87) 


and the fact that the closest ray on the cone to the initial 
desired scatter direction d should^ be chosen. In other words, 

3 0 should be chosen to maximize d-d 0 with the constraint that 
3 0 • e m = cos go* This implies that the following equations need 
to be solved 


3( d-d 0 + xd 0 -e m ) 


= 0 


3 (d.do + Ad Q .e m ) 


= 0 


9 
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where 


d = d x x + dyy + d z z 

A A A A 

d„ = cos<j> sine x + sin* sine v + cose z 
ooo r o o o 


e m = e x x + V + V ' 
This can be solved to give 


cos<{. 0 = 


d x + xe x 


± |( d x + ^x ) 2 +(d y +Ae y ) 2 


sin<j> 0 = 


d + xe 

-1 lx 


± i d x +Xe x ) 2 +(d y +Xe y )2 


cose Q = 


d z + Ae z 


± i d x +Ae x ) 2 +( V Ae y ) 2 +(d z +xe z> 


and 


sine 0 = 


x/(V Ae x> 2 * ( d v * Ae / 
i< d x +Ae x ) 2 + <V xe y> 2 +(d z tAe z> : 


The parameter a is found by substituting these equations into 
Equation (87), which gives 
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A= -(d-e m ) ± (r*e m ) 
m m 


1 - ( d • e ) ‘ 
m 

l-(r-e )‘ 
m 


The signs are chosen such that d-d 0 is maximized. This d 0 and its 
corresponding reflection point and diffraction point (corner point) 
are used as the starting parameters in the incremental method. 

Now that the reflection point and diffraction point are known 
the field can be found. The reflected field off the elliptic 
cylinder is in the form of Equation (2), where the distance from 
the cylinder to the edge is finite. The diffracted field from the 
edge is given by Equation (11), where the incident field on the 
edge is the reflected field off the cylinder, and therefore, the 
wavefront incident on the edge has two radii of curvature. This 
field is given by 


E^d 

cyl ,j ,m 


E 1 (Qn) 


r r 
p l p 2 


(p^+s 1 ) (p 2 +s ' ) 




-jks'e 


-jks 


where s is in the far- field and therefore 


e -jks 

s 


% 


jkR.-d e- JI<d 
e 0 d 


The parameters needed to find the reflected field off the elliptic 
cylinder are given in Equations (70a, b) and (71a-k) and the parameters 
for the diffracted field associated with the edge are given by 
Equations (45a-g) except for the distance parameter L and the 
caustic distances. The incident field on the edge has two radii 
of curvature as stated earlier such that the distance parameter 
given in Equation (14) must be used where s-*» or in other words 

p 1 p 2 sin e 0 


The incident radii of curvature at the edge are due to the radii 
of curvature at the surface of the cylinder plus the distance 
from the point of reflection to the diffraction point or 
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P 1 


=pj + s- 


p 2 


r , 
P o' S 


The caustic pi is the radius of curvature of the incident wave at 
the edge in tne plane of the edge and incident ray and is given 
by Eulers Theorem [66], as 

112 12 

— = — cos a + — cos B . 

p e P 1 p 2 

The angles are defined by the edge and the principal directions of 
the reflected wavefront such that 


cos a = X-, • e„, 
1 m 


cos B = X^ • e . 

2 m 

The principal directions of the reflected wavefront are found from 
Equations (7-9) given in Chapter II. 

The total field for this scattering mechanism is the sum of all 
the contributions from the edges, if they exist, and is given by 


E cyl-pl ^s’^ 


P M. 

I l 

j=l m=l 


^ rd 1 . (0 ,<(>). 

cyl ,j ,m v s y s' 


(88) 


The scattering mechanism that corresponds to the reflected 
field off the cylinder after being diffracted by an edge is il- 
lustrated in Figure 57. The first step in finding the f^eld is 
to determine its ray path for a given scatter direction d. In 
this case the incident vector is dependent on the edge diffraction 
point which in turn is dependent on the reflection point on the 
cylinder. Again, the incremental method must be employed to 
solve this non-linear problem. 
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Figure 57 — II 1 ustration of reflected field off the cylinder after 
diffraction by an edge. 
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The laws of reflection for the reflection off the cylinder 
require that the following equations be satisfied: 

f < v - z 'V* s > = ( "xV n y y (n xy r v d x ) 

+( "xV n y d y )(n x I y- n y I x ) = ° 

and 

g(v,z,9 s ,* s ) = (n x I x +n y I y )d 2+ (n x d x +n y d y )I 2 =0. 

The functions f and g can be expanded in a first order Taylor 
series and solved for 6v and 6z as was done previously in Equations 
(86a and b). The answer is the same and does not need to be 
repeated here. The partial derivatives in this case are given by 


3f , f 9n al an dl \ 

— -1 =( — — - I +n + * I + n (n d -n d ) 

3v l 3v x x 3v 3v y y 3v J ' x y y x' 


/ 3n x 8n v 

+(n x I x +n y I y ) ( 37" d y " a/ d > 


/3n 3 1 3n 

+ (n d +n d )( - — I + n — ^ ^ I 

v v v \/ \//| \/ v 3 \i > 


31. 


- n. 


xx "yy\3v y "x 3v 3v x "y 3v 


an x an v 

. .. . . . — - d + — £ d 
xy yx'l3v x 3v y 


+ (n v I w -n„I ) 


3f j r 8I y 3l v \ r aI v 9 V 

3z " ^ n x d y~ n y d x^ ( n x 3z + n y 3 z~J + ^ n x d x +n y d y^ (°x 3z " n y 3z , 


3f . 


3d. 


3d. 


(n v I +n„I w ) n v — L - n„ )+ (n I -n, I) 


3<j> s l "x‘x "y*y\"x 3<(> s "y 3 <( > S J "xy y x^ x 3<j> s y 3<|> s ^ 


e> 


3d. 


+ n. 
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af. 


30 


= 0 (since sine dependence can be removed from this 
equation) 


and 



3n 31 3n 

x I +n -2L + -X-1 

X X 3V 3V 




+ fn d +n 


31 


v 



n 


5 

x 3z 



KWV 




x 34 


+ n 


3d \ 

—l) i 

y 3^ s J z 



- ("xWJ 




+ n y 



The definitions of n and d and their derivatives have_been given 
before, however, the determination of I, 3l/3v, and 3 1/ 3z need to 
be derived. 

The incident vector (for the ray incident on the cylinder) is 
given by 

T(v,z) = R r (v,z) -R d (v,z) (89) 

where R r is the reflection point and the edge diffraction point is 
given by 
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R d (v.z) - R p + s p 


I(v,z) • e, 


m 


Vl-[I(v,z)*eJ 


e m’ 


(90) 


R_ = C_ + [ ( R_ -C ) • e_]e_ 
p m LX s m' mm 


s = |R - R | . 
p 1 s p 1 


This leads to the equation 


I = (R -RJ - j==L e 


r p' 






m 


(91a) 


and 


1 1 1 2 = l* I = l R r - R p i 


■ 2 + 2 lie) 


J 


:x2 


l-(I-e) 


P l-(I-e) 2 


(91b) 


Taking the dot product of Equation (91a) and e m and then squaring 
the equation gives 


iv, 2 


(1-0 » ~2s [(R -R ) *e ] 


(I-eJ 


m' 


m 


r p' nr p 1 -' r p' m- 


I /s A 


slu-if 

p' nr 

A /X 

1 " ( 1 " 6 ) 

' m 


2 . 


(92) 
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Now, substituting Equation (91b) into Equation (92) then squaring 
the equation again to remove the square root term, the equation 
can be regrouped in the form 

0 o a 4 + C-,cx 2 + C = 0 
2 1 o 


where 


C 2 < s p-IV R p| Z ) 2 + 4s p[<W 4 / 
C 1 - - 2 < s p + IVV 2 > 

c 0 ■ r<W 


This can be solved using the quadratic formula solution and is 
given by 


(I 


e m ) 


= + 


- C 1 ± 


C 2 - 
L 1 


4C o C 2 


2 C, 


(93) 


This can be used in Equations (90) to find the diffraction point 
and hence in Equation (89) to find the incident vector. The signs 
are chosen in Equation (93) to minimize the distance from the source 
to the edge to the reflection point. 

The derivatives of the incident vector are found next. The 
normalized incident vector is given by 


m 

so 


160 




Taking the dot product with respect to the edge vector e m of 
Equation (96) and then rearranging gives 



Substituting Equation (97) into Equation (95) gives the solution 
for the partial derivatives of the incident vector. These can be 
used to find 6v and 6z and hence the reflection point and diffraction 
point using the same prodecure discussed for the previous scattering 
mechanism. 

A known starting location for the first computation is found 
next. This can be found by defining a fixed point on the cylinder 
and from this define the diffraction point and scatter direction 
d 0 for this point. The best location to define a known point R 0 
on the cylinder for this problem is a point on the end cap rim. 

The diffraction point is given by 
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where 


R . = R n + t c e m 
d p s m 


t = (R - RJ • e 
e o s m 

s s - IT* - V 

s o " l», - V' 

/x A 

R„ = [(R c -C m ) *e m ]e m + C m 
p L ' s m' m J m m 

A A 

R„ rt = [(R„ - C m ) • el e m + C. 
po o m m m m 


This is illustrated in Figure 58. When the diffraction point is 
found then the scatter direction can be found from the following 
equations: 


I = 


R o -_ R d 

IV R dl 


<f> = tan 

so 


-1 


(n 2 x -n^)I y - 2 n A l x 


9 


SO 


cos" 1 (I z ) 
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d 0 = cos<j> S g sin e so x + sin$soSinesgy+coses 0 z- This d o and its 
corresponding reflection point (R 0 ) and diffraction point are used 
as the starting parameters in the incremental method. 



Figure 58 — 111 ustration of starting point for finding diffraction 
point on an edge and the reflection point on the 
cylinder for the diffraction-reflection problem. 


Once the reflection point and diffraction points are determined 
the field can be found. The diffracted field of the edge is in the 
form of Equation (11), where the distance from the edge to the 
cylinder is finite, the reflected field of the cylinder is given 
by Equation (2), where the incident field on the cylinder is the 
diffracted field from the edge. This incident field has two 
principal radii of curvature. The diffracted-reflected field is 
given by 
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E dr 

j ,m,cyl 


= E 


(Q e ) 


R s/s"(s'+s") 


J 


r r 
p l p 2 


-jks' 


-jks 


where s is in the far field and therefore 

e -jks ^ jk R r -d e -jkd 
i e d . 

The parameters needed to find the diffracted field from an edge 
are given in Equations (45a-g), except the distance parameter needed 
here is in the near field and is given by 


where s' is the distance from the source to the diffraction point 
and s" is the distance from the diffraction point to the reflection 
point. The parameters for the reflected field are given in 
Equations (71a-k) except the incident field has two radii of curvature 
so the reflection caustic distances need to be found from Equation 
(5) with which is given by 



where 

i A 

p-j = s , in the principal direction <j> 

p 2 = s"+s', in the principal direction B . 

The parameter cose 1 is found in Equation (71c), R, in Equation 
(71b), and e.., |e|in Equations (6a and b) . 

* J 
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The total field for this scattering mechanism is the sum of 
all the contributions from the edges if they exist and is given by 


P M. 

E pl -cyl ^ e s 5< M “ 2 ^ m | 1 E j,m,cyl^ e s**s^ 


(98) 


The total field for the aircraft model is the sum of all the 
scattered fields^ from the finite plates, listed in Table II, 
represented by E p -| ; from the finite elliptic cylinder, listed in 
Table III, represented by E^yi ; and from the plate-cylinder inter- 
action terms, listed in Table IV. This total field can therefore 
be written as 


E(0 »<(>) - E^ + Ei + E -1+ I" i + E^'? 

' s v s pi cyl pi -cyl pi -cyl 


+ F^ r + + E^ r 

t cyl-pl L cyl -pi L pl-cyl * 


TABLE IV 

SCATTERED FIELDS INCLUDED INVOLVING PLATE-CYLINDER INTERACTIONS 



Symbol 

Description 

Figure 

Equation 

1) 

E rr 

pl-cly 

Reflected field off the 
cylinder after reflection 
by a plate 

52 

81 

2) 

fTt 
pi -cyl 

Transition field off the 
cylinder after reflection 
by a plate 

53 

82 

3) 

-Err 

t cyl-pl 

Reflected field off a 
plate after reflection 
by the cylinder 

54 

83 

4) 

E rd 

L cyl-pl 

Diffracted field off an 
edge after reflection 
by the cylinder 

55 

88 

5) 

E dr 
pi -cyl 

Reflected field off the 
cylinder after diffraction 
by an edge 

57 

98 
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This result, then, includes the principal contributors to the 
scattered fields from the finite plates, finite elliptic cylinder, 
and plate-cylinder interaction terms. They have been found by 
first tracing the ray paths for a given scatter direction, then 
testing these ray paths to see whether or not the fields are 
shadowed by any obstacle along the ray path. The non-shadowed 
fields are then computed using geometrical optics and/or the 
geometrical theory of diffraction and superimposed to give the 
total field. 

To illustrate the effects that some of these scattering mech- 
anisms have on the radiation pattern of an antenna mounted near a 
plate-cylinder structure, a simple model consisting of a flat 
plate attached to a finite elliptic cylinder in the y-z plane, 
as shown in Figure 59, is studied. The elevation plane (x-z 
plane) radiation pattern for an infinitesimal slot mounted parallel 
to the y-axis is shown in Figure 60a. The cylinder has little 
effect on the radiation pattern in this plane. The roll plane 
(x-y plane) radiation pattern for the same source is shown in 
Figure 60b. The need for a corner diffraction term is shown by 
the discontinuity in the right side of the backlobe of the pattern. 

The result using the recently developed corner diffraction term is 
shown as the dotted line in the same figure. The discontinuity in 
the left side of the backlobe is due to shadowing by the cylinder. 

The mechanism to compensate for this discontinuity has not been 
studied at this time, but it will have only a small effect other 
than the discontinuity. The roll plane radiation pattern for an 
infinitesimal slot located parallel to the z-axis is shown in 
Figure 60c. The lack of a reflected-diffracted term is apparent 
by the large discontinuity in the solid line at <j> = 90°. The 
smoothing out of the discontinuity by including the reflected- 
diffracted term is shown by the dotted line. Even for this 
simple structure the interactions between the different parts are 
apparent in their area of greatest influence. The way that the 
GTD compensates for these discontinuities, however, is also quite 
apparent. This is useful to keep in mind as the structures become 
more complicated, and consequently, the interaction terms become 
more involved and less practical to include in the solution. 

To test the validity of the approximations involved in the finite 
elliptic cylinder aircraft model, an experimental model that 
resembles the computer model was built and tested. The model 
consisted of a finite cylinder with two flat plates resembling 
wings as shown in the insert of Figure 61. The comparison of 
calculated and measured results for the roll plane radiation 
pattern is shown in Figure 61, The overall agreement between the 
two results is very good, even down to some of the small back" 
lobes. The disagreement at <j>= -90° is due to the lack of a 
diffracted-reflected-diffracted term. This third-order term 
is hard to calculate and seems to be important only in this 
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Figure 59--Geometry of a preliminary model. 
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Figure 60a--Elevation plane radiation pattern ( E© ) of small slot 
mounted parallel to y-axis on plate. 
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Figure 60b--Roll plane radiation pattern (Eg) of small slot 
mounted parallel to y-axis on plate, showing 
effect of corner diffraction term. 
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Figure 60c— Roll plane radiation pattern (E^) of small slot 
mounted parallel to z-axis on plate, showing 
effect of reflected-diffracted term. 
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particular pattern cut. Therefore, it does not seem necessary to 
include this term at this time. The discontinuity at <j> = 120° is 
due to some higher order terms involving the second wing which 
is not included. Again, these do not appear to be of engineering 
significance so they have been left out of the solution. This 
comparison indicates that the theory can predict patterns for 
simple structures resembling the analytic model. 

It is interesting, at this point, to compare the results in 
Figure 61 with the flat plate computer representation of the 
simple aircraft model as shown in the insert to Figure 62. This 
emphasizes the difficulties involved in representing a smooth 
convex surface by a flat plate as was previously used for the F-4 
aircraft in Figures 33a-f. The radiation pattern for the flat 
plate version of the simple aircraft model is shown in Figure 62. 

The major features of the patterns are quite different. This is 
simply due to the fact that the specular reflection off the flat 
plate does not adequately approximate the specular reflection off 
of the cylinder as might be expected. This flat plate model can 
then be expected only to give a crude approximation to an aircraft 
only if the curved surfaces are not strongly illuminated. 

The next question that needs to be resolved is how well does the 
finite elliptic cylinder model predict the radiation patterns of 
antennas mounted on real aircraft shapes. The finite elliptic 
cylinder aircraft model has been used to compute some of the 
patterns that were measured at NWC using a single slot on the wing 
of the F-4 scale model as discussed in Chapter III. The geometry 
of this theoretical model is shown in Figure 63. The dimensions 
for the model are listed in Table V. The antenna locations are 
shown in Figure 32. A comparison of the measured and calculated 
radiation patterns for the main beam elevation plane for the slot 
at location #2 (refer to Figure 31a) is shown in Figure 64a. This 
result can be compared with the flat plate result that was used to 
crudely model the aircraft in Figure 33a. Note that the ripple 
level, in the region <j> = 60° to 150° , due to the reflections off 
the fuselage has been improved as expected. The same problem 
of comparing near field measured results and far field computed 
results still exists in all these patterns. However, a few of the 
results for the various cases will be shown here for illustrative 
purposes. In each case, Figures 31 a-d show the aircraft pattern 
orientation. The polarization of the received signal is computed 
in the same way the measurements were taken, that is, with the 
probe aligned with the antenna to give the dominant polarization. 
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Figure 61 --Comparison of measured and calculated results for 
the roll plane radiation pattern (E<j,) of a simple 
aircraft model with a short monopole mounted on 
its wing, (freq = 10.87 GHz). 



* 



Figure 62--Roll plane radiation pattern (E^) of a flat plate 
representation of the simple aircraft model with 
a short monopole mounted on its wing. 
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(a) FRONT VIEW 
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Figure 63 — Illustration of geometry of F-4 aircraft model used 
in finite elliptic cylinder model. 



TABLE V 

DIMENSIONS FOR FINITE ELLIPTIC CYLINDER MODEL OF F-4 
All Dimensions in Inches for Full Scale Model 


Fuselage 


fl Top ' 87 • 
A Bottom = 19 ‘ 


B = 54. 


Z End Cap Nose " -420 *’ a Nose ~ 90 ° 
Z End Cap Tail = 26 - 6 ’ a Tail = 20 ° 


Wings 
Plate 1 


Corner 

X 

Y 

Z 

1 

0 . 

54. 

-228. 

2 

0 . 

166.2 

- 86.2 

3 

0 . 

166.2 

32.8 

4 

0 . 

54. 

0 . 



Plate 2 



Corner 

X 

Y 

Z 

1 

0 . 

166.2 

- 86.2 

2 

0 . 

166.2 

- 96.2 

3 

19.7 

228.7 

0 . 

4 

19.7 

228.7 

53.4 

5 

0 . 

166.2 

32.8 
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TABLE V cont. 



Plate 3 


Corner 

X 

Y 

1 

0 . 

-54. 

2 

0 . 

-166.2 

3 

0 . 

-166.2 

4 

0 . 

-54. 


Plate 4 


Corner 

X 

Y 

1 

0 . 

-166.2 

2 

19.7 

-228.7 

3 

19.7 

-228.7 

4 

0 . 

-166.2 

5 

0 . 

-166.2 


Plate 5 


Corner 

X 

Y 

1 

87. 

0 . 

2 

149.9 

0 . 

3 

154.3 

0 . 

4 

87. 

0 . 


Source Coordinates 

Source 

X 

Y 

1 

5.14 

182.5 

2 

0.00001 

149.1 

3 

0.00001 

126.3 

4 

0.00001 

150.3 

5 

0.00001 

111. 


Z 

0 . 

32.8 

- 86.2 

-228. 


Z 

32.8 

53.4 

0 . 

-96.2 

- 86.2 


Z 

205.1 
197.3 

153.1 
11.4 


Z 

-36.91 

-73.6 

-102.5 

-4.6 

-74.1 




The comparison of the main beam elevation plane pattern with 
a slot antenna at location #4 is shown in Figure 64b. This result 
can be compared with Figure 33b. The pattern in the main beam roll 
plane (Figure 31c) for a slot antenna at location #2 is shown in 
Figure 64c and compares with Figure 33c. The elevation plane 
pattern (Figure 31b) for a slot antenna at location # 2 is shown in 
Figure 64d and compares with Figure 33d. Its roll plane pattern 
(Figure 31a) is shown in Figure 64e and compares with Figure 33e. 

The roll plane result especially emphasizes the specular reflection 
off the fuselages curved surface. The elevation plane result for 
a slot antenna oriented such that its maximum radiation points 
aft of. the aircraft at location #5 is shown in Figure 64f and 
compares with Figure 33f. 

In all of the above results, the major characteristics of the 
patterns do show some agreement. In this model the horizontal 
stabilizers have little effect so they are neglected for this problem. 
Some higher interaction terms that were not included are significant 
in small angular regions in some of the various pattern cuts; how- 
ever, it is not deemed feasible to include them at this time. 

In order to further test the finite elliptic cylinder model's 
ability to predict radiation patterns of antennas mounted on the 
wings of real aircraft, pattern measurements of a slot antenna 
mounted on the wing of a 1/20 scale model of a Boeing 737 were taken 
at NASA (Hampton, Va.). The geometry of the finite elliptic 
cylinder model of the Boeing 737 is shown in Figure 65. The 
dimensions used for this model are given in Table VI with a fre- 
quency of 1.75 GHz, which corresponds to 35 GHz for the scale model. 
The source is a Ka band waveguide mounted parallel to the y-axis 
in Figure 65 with the location given in Table VI. The patterns 
are presented with respect to the coordinate system in Figure 66. 

The Eq> radiation patterns are shown in Figures 67a-e for patterns 
cut at angles <p' = 0° to <j> 1 = 150° in steps of 30°, respectively. 

The patterns have been normalized to the maximum of the pattern 
in the <j>' = 0° cut. The pattern at tj>' = 90° has not been shown 
since for this pattern cut the amplitude of the E 0 ' polarized 
field is below -30 dB which is in the noise region for the meas- 
ured pattern. In all of the patterns shown, the finite elliptic 
cylinder model tends to predict very well the frequency and level 
of the ripple in the patterns in most regions. Some of the dis- 
crepancies can be attributed to the following reasons. The high 
frequency of the ripple in the patterns means that any misalignment 
in the experimental measurements could cause shifts in the actual 
placement of the peaks and nulls. Also, the support structure 
appears to cause some scattering in the e 1 = 150° to 210° region 
in all the patterns. The measurements in the <p‘ = 30° and <j>' = 150° 
pattern cuts in the e' = -90° and e 1 = 90° regions have pointed 
out that the nose and tail configurations of the finite elliptic 
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cylinder model are not completely adequate to predict all the 
detail accurately. However, the discrepencies do not appear to be 
severe. 

Overall, the agreement between the calculated and (NWC and 
NASA) measured results appears to be good for engineering appli- 
cations. They both point out that the fuselage can strongly 
effect the patterns of antennas mounted on the wings of aircraft. 

In particular the reflections off the fuselage and other higher 
order interactions tend to cause a lot of ripple in the patterns. 

This can be detrimental to system performance in many instances. 
However, if the antennas are made directive such that they do not 
illuminate the fuselage strongly and the antenna locations are chosen 
judiciously the ill effects of the fuselage and other scatters can 
be minimized. Thus using this numerical solution, the antenna 
designer can begin to consider antenna locations on the aircraft 
wing without having to make numerous scale model measurements. 

As a practical example of how the finite elliptic cylinder 
model can be used to solve an aircraft antenna design problem, the 
design of an antenna system to give coverage in the aft reqion of 
an F-4 aircraft for quadrant detection at UHF, as undertaken in 
Reference [33], is briefly discussed. Classical array analysis was 
used to find the element spacing, scan angle, and distribution of 
a two-element array to obtain a desired pattern shape above the 
wing. The array was then located on the wing using the finite 
elliptic cylinder model to optimize the antenna pattern. The optimum 
location was found to be at location #2 of Figure 32 with the array 
being mounted at 45° with respect to the fuselage. The volumetric 
pattern predicted for this antenna system on the F-4 is plotted in 
Figure 68 in terms of its directive gain with respect to the co- 
ordinate system in Figure 66. The desired region of coverage is shown 
as the rectangular box in the figure. The directivity, which is the 
maximum directive gain, for the antenna system in place on the 
aircraft was determined to be 12.3 dB. The directive gain is 
plotted in Figure 68 using different colors to represent the 
different gain levels. Each change in color represents a 5 dB window 
showing regions of farther reduction in the gain. The volumetric 
pattern shows that the two-element array design does a realistic 
job of approximating the desired objectives. Other optimum 
antenna designs can be accomplished in a similar efficient 
manner on a digital computer using the finite elliptic cylinder 
model and then tested experimentally to verify the findings. 
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Figure 64a--Comparison of NWC measured results on an F-4 model 
with the finite elliptic cylinder model results of 
main beam elevation plane radiation pattern for 
antenna location #2. 


179 





Figure 64b — Comparison of measured and calculated results of 
main beam elevation plane pattern for antenna 
location #4. 
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Figure 64d— Comparison of measured and calculated results of 
elevation plane pattern for antenna location #2. 
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Figure 64e--Comparison of measured and calculated results of 
roll plane pattern for antenna location #2. 
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Figure 64f--Comparison of measured and calculated results of 
of elevation plane pattern for antenna location 
#5 with the slot pointed straight aft of aircraft. 
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TABLE VI 

DIMENSIONS FOR FINITE ELLIPTIC CYLINDER MODEL OF BOEING 737 
All Dimensions in Inches for Full Scale Model 
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\y* 

Figure 66--Coordinate system used for volumetric patterns. 
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Figure 67a--Comparison of NASA's measured E 0 i results on a 

Boeing 737 model with the finite elliptic cylinder 
model results at <t>' =0. 
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Figure 67c--Comparison of measured and calculated Ee' results 
. for <j>' = 60°. 
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Figure 67d— Comparison of measured and calculated Ee> results 
for = 120°. 
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Figure 68. Volumetric directive gain results with various gain levels represented in different 
colors for optimum two element array design mounted on the wing of an F-4'model, 



CHAPTER VI 

SUMMARY AND CONCLUSIONS 


The object of this research has been to analyze radiation 
patterns for aircraft wing mounted antennas. This is an 
analytical study of basic antenna types for frequencies above about 
100 megahertz in which the antenna is mounted on or near a wing. 

The aircraft is modelled in its most basic form so that the study 
is applicable to general -type aircraft. The fuselage is modelled as 
a perfectly conducting finite elliptic cylinder. The wings and 
horizontal and vertical stabilizers are modelled as perfectly con- 
ducting "n" sided flat plates that can be arbitrarily attached to the 
fuselage or to themselves. Presently, the antenna locations are 
assumed to be on the upper surfaces of the wings or at locations 
removed from engines and stores such that these effects are negligible. 
The theory, however, is applicable to the addition of engine and 
store effects by modelling these structures as flat plates or small 
finite elliptic cylinders. 

The basic approach applied in this study has been to separately 
analyze the basic scattering shapes making up the computer simu- 
lation model of an aircraft. The single structures are analyzed 
using the ray optics approach such that only the dominant effects 
are included. The scattered fields of the various component parts 
of the structure are found using the Geometrical Theory of Diffraction. 
The basic solutions needed for this study were briefly discussed in 
Chapter II. This includes the geometrical optics reflected fields 
from curved surfaces, the diffraction by a wedge and curved wedge, 
and the slope diffraction by a wedge using arbitrarily oriented 
dipole moment sources. A heuristic corner diffraction coefficient 
was, also, introduced. Curved surface diffraction was presented 
along with a newly developed transition function which describe 
the scattered field behavior around the shadow boundaries of 
curved surfaces. 

The scattering from a number of finite flat was discussed 
in Chapter III. The plates can have an arbitrary orientation 
in space and can be attached to form bent plates. The present 
solution is restricted in that only first and second order terms 
are included except for double diffraction which is normally neg- 
ligible. This solution is useful for antennas mounted near a flat 
plate, corner reflectors, and even a crude flat plate aircraft. 

These plates are used to model wings and horizontal and vertical 
stabilizers on the finite elliptic cylinder model of an aircraft. 

They can, also, be used to model other scattering objects such 
as buildings or ships with some modification. 
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The scattering from a finite elliptic cylinder was discussed 
in Chapter IV. The end caps of the cylinder can form an angle with 
the cylinder axis. The scattered fields are found using geomerical 
optics, the Geometrical Theory of Diffraction and a curved surface 
transition function. The validity of the solution was shown by 
comparing the results against exact and moment method solutions 
and experimental results. The finite elliptic cylinder is used to 
model the fuselage of an aircraft and can be used to model engines 
and stores. The solution can also be used to model other objects 
such as ship masts or hulls. 

The combination of the above scattering objects to form a 
general -type aircraft model is discussed in Chapter V. Measured 
results for a simple aircraft-like model that resembles the finite 
elliptic cylinder model verified the approximations in the 
theoretical solutions. Measurements taken at NWC and NASA on 
actual aircraft models verified the usefulness of the finite 
elliptic cylinder model in predicting the radiation patterns for 
wing-mounted antennas. 

In general, the results show that an antenna mounted on a 
wing can interact strongly with the fuselage and wings causing 
large amounts of ripple in the pattern. This suggests that proper 
element design especially directive antennas that do not illuminate 
the fuselage strongly, and proper placement of the antenna on the 
wings. to minimize the interference problems are important con- 
siderations that must be taken into account when designing wing 
mounted antenna systems. The finite elliptic cylinder model 
should prove to be a very useful tool to an antenna designer in 
anticipating these problems in a very economical and efficient 
manner. The computer code used to generate a pattern for one 
antenna element typically takes only 30 seconds or less on 
an IBM 370 digital computer. 

A brief discussion of a practical antenna design problem is 
presented to show the usefulness of the finite elliptic cylinder 
model. The model was used to optimize the design, of a quadrant 
detection antenna mounted on the wing of an F-4 aircraft. A 
volumetric pattern of the predicted antenna systems performances 
illustrated how closely the antenna design can come to the ideal 
characteristics of a quadrant detection antenna. 

The versatility of the present solution can be increased 
in many ways. The engine and store effects can be added to the 
aircraft model by using small finite elliptic cylinders or plates 
to represent these objects. The comparisons of model measurements 
with the finite elliptic cylinder have indicated that the nose 
and tail representation are not always adequate in some small 
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angular regions. This could be improved by including equivalent 
currents to remove caustic effects and perhaps by alternative 
representations of the nose and tail. The solution can also be 
extended to allow the analysis of tail mounted antennas. 


196 



REFERENCES 


[1] Carter, P. S., "Antenna Arrays Around Cylinders," Proc. IRE, 
Vol. 31, December 1943, pp. 671-693. 

[2] Carter, P. S., "Antenna and Cylindrical Fuselage," Report 
No. 895-11, RCA Laboratories, Rocky Point, N.Y.; December 24, 
1943. 

[3] Sinclair, G., "The Patterns of Antennas Located Near Cylinders 
of Elliptical Cross Section," Proc. IRE, Vol. 39, No. 6, 

dune 1951 . 

[4] Richmond, J. H., "A Wire-Grid Model for Scattering by Con- 
ducting Bodies," IEEE Trans. Antennas and Propagation, 

Vol. AP-14, pp. 782-786, November 1966. 

[5] Richards, G. A., "A Boundary-Value Technique for Computing 
the Patterns of an Antenna Near a Conducting Body of Arbitrary 
Shape," Report 2235-2, 11 September 19767, The Ohio State 
University ElectroScience Laboratory, Department of Electrical 
Engineering; prepared under Contract N10019-67-C-0063 for 
Naval Air Systems Command. 

[6] Richmond, J. R. , "Computer Analysis of Three-Dimensional Wire 
Antennas," Report 2708-4, 22 December 1969, The Ohio State 
University ElectroScience Laboratory, Department of Electrical 
Engineering; prepared under Contract DAAD 05-69-C-0031 for 
Ballistic Research Laboratory. 

[7] Wang, N. N., Richmond, J. H., and Gilreath, M. C., "Sinusoidal 
Reaction Formulation for Radiation and Scattering from Con- 
ducting Surfaces," IEEE Trans, on Antennas and Propagation, 
Vol. AP-23, No. 3 (May 1975), pp. 376-382. 

[8] Ryan, C. E., Jr., and Rudduck, R. C., "Calculation of the 
Radiation Pattern of a General Parallel -PI ate Waveguide 
Aperture for the TEM and TEg-, Waveguide Modes," Report 1693-4, 
10 September 1964, The Ohio State University ElectroScience 
Laboratory, Department of Electrical Engineering; prepared 
under Contract N62269-2183 for U. S. Naval Air Development 
Center, Johnsville, Pennsylvania. 

[9] Ryan, C. E., Jr., and Rudduck, R. C., "A Wedge Diffraction 
Analysis of the Radiation Patterns of Parallel-Plate Wave- 
guides," IEEE Trans, on Antennas and Propagation Comm., 

Vol. AP-16, No. 4, July 1968. 


197 



[10] Wu, D. C. F., Rudduck, R.C., and Pelton, E. L., "Application 
of a Surface Integration Technique to Parallel-Plate Wave- 
guide Radiation Pattern Analysis," IEEE Trans. Antennas and 
Propagation, Vol . AP-17 (May 1969), pp. 280-285. 

[11] Burnside, W. D. and Pelton, E. L. , "Wedge Diffraction Theory 
Analysis of Parallel-Plate Waveguide Arrays," Report 2382-14, 
25 September 1969, The Ohio State University ElectroScience 
Laboratory, Department of Electrical Engineering; prepared 
under Contract F336 1 5-67-C-l 507 for the Air Force Avionics 
Laboratory at Wright-Patterson Air Force Base, Ohio. 
(AFAL-TR-69-240) (AD 860082) 

[12] Burnside, W. D. , Pelton, E. L. , and Peters, L., Jr., 

"Analysis of Finite Parallel-Plate Waveguide Arrays," to be 
published as a Communication in IEEE Trans, on Antennas and 
Propagation. 

[13] Yu, J. S., Rudduck, R. C. , and Peters, L., Jr., "Comprehensive 
Analysis for E-Plane of Horn Antennas for Edge Diffraction 
Theory," IEEE Trans, on Antennas and Propagation, Vol. AP-14, 
(March 1966) pp. 138-149. 

[14] Yu, J. S., and Rudduck, R. C., "H-Plane Pattern of a Pyramidal 
Horn," IEEE Trans, on Antennas and Propagation Comm., Vol. 
AP-17, No. 5, September 1969. 

[15] Ratnasiri, P. A. J., Kouyoumjian, R. G. , and Pathak, P. H., 
"The Wide Angle Side Lobes of Reflector Antennas," Report 
2183-1, 23 March 1970, The Ohio State University Electro- 
Science Laboratory, Department of Electrical Engineering; 
prepared under Contract AF 1 9 ( 628) -5929 for Air Force 
Cambridge Research Laboratories. 

[16] Mentzer, C. A., Pathak, P. H., and Peters, L., Jr., "Pattern 
Analysis of an Offset Fed Parabolic Reflector Antennas," 

Report 3220-2, June 1972, The Ohio State University Electro- 
Science Laboratory, Department of Electrical Engineering; 
prepared under Contract N001 78-71 -C-0264 for U. S. Naval 
Weapons Laboratory. 

[17] Ryan, C. E., Jr., and Rudduck, R. C., "Radiation Patterns 
of Rectangular Waveguides," IEEE Trans, on Antennas and 
Propagation Comm., Vol. AP-16 (July 1968), pp. 488-489. 

[18] Balanis, C. A., and Peters, L., Jr., "Analysis of Aperture 
Radiation from an Axially Slotted Circular Conducting Cylinder 
Using Geometrical Theory of Diffraction," IEEE T-AP-17, 

January 1969, pp. 93097. 


198 



[19] Balanis, C. A., and Peters, L., Jr., "Aperture Radiation from 
an Axially Slotted Elliptical Conducting Cylinder Using 
Geometrical Theory of Diffraction," IEEE T-AP-17, July 1969, 
pp. 507-513. 

[20] Balanis, C. A., and Peters, L., Jr., "Radiation from TE-|q 
M ode Slots on Circular and Elliptical Cylinders," IEEE 
T-AP-18, No. 3, May 1970, pp. 400-403. 

[21] Ryan, C. E., Jr., "Analysis of Radiation Patterns of Antennas 
on Finite Circular Cylinders and Conically-Capped Cylinders," 
Report 2805-2, 25 September 1970, The Ohio State University 
ElectroScience Laboratory, Department of Electrical Engineer- 
ing; preapred under Contract DAA21-69-C-0535 for Picatinny 
Arsenal. (AD 875581L) 

[22] Burnside, W. D., Gilreath, M. C., Marhefka, R. J., and 
Yu, C. L., "A Study of KC-135 Aircraft Antenna Patterns," 

IEEE Trans, on Antennas and Propagation, Vol . AP-23, No. 3, 

May 1975, pp. 309-16. 

[23] Burnside, W. D. , "Analysis of On-Aircraft Antenna Patterns," 
Report 3390-1, August 1972, The Ohio State University Electro- 
Science Laboratory, Department of Electrical Engineering; 
prepared under Contract N62269-72-C-0354 for Naval Air 
Development Center. (AD 777 989) 

[24] Yu, C. L., and Burnside, W. D. , "Volumetric Pattern Analysis 
of Fuselage Mounted Airborne Antennas," Report 2902-24, 

April 1976, The Ohio State University ElectroScience 
Laboratory, Department of Electrical Engineering; prepared 
under Grant NGL 36-003-138 for National Aeronautics and Space 
Administration, Langley Research Center. 

[25] Marhefka, R. J., "Roll Plane Analysis of On-Aircraft Antennas," 
Report 3188-1, December 1971, The Ohio State University 
ElectroScience Laboratory, Department of Electrical Engineer- 
ing; prepared under Contract N62269-71 -C-0296 for Naval 

Air Development Center. 

[26] Yu, C. L., and Burnside, W. D. , "Elevation Plane Analysis 
of On-Aircraft Antennas," Report 3188-2, January 1972, The 
Ohio State University ElectroScience Laboratory, Department 
of Electrical Engineering; prepared under Contract N62269- 
71-C-0296 for Naval Air Development Center. 

[27] Burnside, W. D., Marhefka, R. J., and Yu, C. L. , "Roll Plane 
Analysis of On-Aircraft Antennas," IEEE Trans, on Antennas 
and Propagation, Vol. AP-21, November 1973, pp. 780-786. 


199 



[28] Marhefka, R. J., and Burnside, W. D. , "Numerical Solutions 
to Some On-Aircraft Antenna Pattern Problems," Technical 
Report 3390-4, October 1973, Department of Electrical 
Engineering; prepared under Contract N62269-72-C-0354 for 
Naval Air Development Center. (AD 777977) 

[29] Yu, C. L., and Burnside, W. D. , "MLS Airborne Antenna 
Research," Semi-Annual Report 2902-22, May 1975, The Ohio 
State University ElectroScience Laboratory, Department of 
Electrical Engineering; prepared under Grant NGL 36-003-138 
for National Aeronautics and Space Administration, Langley 
Research Center. (NASA-CR-145398) (N 76-10378) 

[30] Yu, C. L. , and Burnside, W. D., "Research on MLS Airborne 
Antenna," Semi-Annual Report 2902-23, April 1976, The Ohio 
State University ElectroScience Laboratory, Department of 
Electrical Engineering; prepared under Grant NGL 36-003-138 
for National Aeronautics and Space Administration, Langley 
Research Center. 

[31] Marhefka, R. J., and Burnside, W. D., "Vertically Polarized 
IFF Antennas Mounted on the Leading Edge of Wings," Report 
2836-2, August 1970, The Ohio State University ElectroScience 
Laboratory, Department of Electrical Engineering; prepared 
under Contract N62269-69-C-0533 for Naval Air Development 
Center. 

[32] Yu, C. L., "Wedge Diffraction Analysis of Antennas of 
Simulated Wing Models," Report 2836-3, September 1970, 

The Ohio State University ElectroScience Laboratory, 
Department of Electrical Engineering; prepared under 
Contract N62269-69-C-0533 for Naval Air Development 
Center. 

[33] Marhefka, R. J. and Burnside, W. D. , "Pattern Optimization 
of Wing-Mounted Array Antennas at UHF," Report 3884-1, 

January 1976, The Ohio State University ElectroScience 
Laboratory, Department of Electrical Engineering; prepared 
under Contract N001 63-74-C-l 580 for Naval Regional 
Procurement Office, Long Beach, California. 

[34] Kouyoumjian, R. G. , and Pathak, P,, "A Uniform Geometrical 
Theory of Diffraction for an Edge of a Perfectly Conducting 
Surface," Proc. of the IEEE, Vol. 62, No. 11, November, 

1974, pp. 1448-1461. 

[35] Deschamps, G. , "Ray Techniques in Electromagnetics," Proc. 
IEEE, Vol. 60, pp. 1022-1035, September 1972. 


200 



[36] Sommerfeld, A., Optics, Academic Press, Inc., New York, 

1954, pp. 245-265] 

[37] Keller, J. B., "Geometrical Theory of Diffraction," Jour. 
Optical Soc. Amer. , _52 (February 1962), pp. 116-130. 

[38] Rudduck, R. C., "Application of Wedge Diffraction to Antenna 
Theory," Report 1691-13, 30 June 1965, The Ohio State Univer- 
sity ElectroScience Laboratory, Department of Electrical 
Engineering; prepared under Grant NSF-338 for National 
Aeronautics and Space Administration, Washington, D. C. 

Also published as NASA Report CR-372. 

[39] Pauli, W., "An Asymptotic Series for Functions in the Theory 
of Diffraction of Light," Phys. Rev., 54 (1 December 1938), 
pp. 924-931. 

[40] Hutchins, D. L., and Kouyoumjian, R. G., "A New Asymptotic 
Solution to the Diffraction by a Wedge," URSI 1967 Spring 
Meeting, Ottawa, Canada, pp. 154-155. 

[41] Hutchins, D. L., "Asymptotic Series Describing the Diffrac- 
tion of a Plane Wave by a Two-Dimensional Wedge of Arbitrary 
Angle," Ph.D. Dissertation, The Ohio State University, 
Department of Electrical Engineering, 1967. 

[42] Pathak, P. H., and Kouyoumjian, R. G., "The Dyadic 
Diffraction Coefficient for a Perfectly Conducting Wedge," 
Scientific Report No. 5, Report 2183-4, 5 June 1970, The 
Ohio State University ElectroScience Laboratory, Department 
of Electrical Engineering; prepared under Contract 
AF19(628)-5929 for Air Force Cambridge Research Laboratories. 
(AFCRL-69-0546) (AD 707821) 

[43] Rudduck, R. C. and Wu, D. C. F., "Slope Diffraction Analysis 
of TEM Parallel-Plate Guide Radiation Patterns," IEEE Trans. 
Antennas Propagat. (Commun.), Vol. AP-17, pp. 797-799, 
November 1969. 

[44] Karp, S. N. and Keller, J. B.: Optica Acta £3, 61 (1961). 

[45] Mentzer, C. A., Peters, L., Jr., and Rudduck, R. C., "Slope 
Diffraction and Its Application to Horns," IEEE Trans. 
Antennas Propagat., Vol. AP-23, No. 2, pp. 153-158, 

March 1975. 

[46] Hwang, Y. M. and Kouyoumjian, R. G. , "A Dyadic Diffraction 
Coefficient for an Electromagnetic Wave Which is Rapidly 
Varying at an Edge," USNC-URSI 1974 Annual Meeting, Boulder, 
Colorado, October 1974. 


201 



[47] Hwang, Y. M. and Kouyoumj i an , R. G., "On the Dyadic Slope 
Diffraction Coefficient," USNC-URSI Meeting, Urbana, 
Illinois, June 1975. 

[48] Walter, C. H., Traveling Wave Antennas , Dover Publications, 
Inc., New York, 1970, pp. 15-16. 

[49] Satterwhite, R. , "Diffraction by a Quarter Plane, the Exact 
Solution, and Some Numerical Results," IEEE Trans. Antenna 
and Propagation, Vol . AP-22, No. 3, pp. 500-502, May 1974. 

[50] Burnside, W. D. , personal communications. 

[51] Keller, J. B., "Diffraction by a Convex Cylinder," Trans. 
IRE, Vol. AP-24, 1956, pp, 312-321. 

[52] Franz, W. , and Depperman, K. , "Theorie der Baugang Am 
Zylinder unter Berucksichtigung der Kreichwel 1 e," Am. 

Physik, Vol. 10, June 1952, pp. 361-373. 

[53] Franz, W., and Depperman, K. , "Theorie der Beugang der 
Kugel unter Berucksichtigung der Kreighwel le," Ann. Physik, 
Vol. 14, June 1954, pp. 253-264. 

[54] Kouyoumjian, R. G. , "Asymptotic High-Frequency Methods," 
Proc. of the IEEE, Vol. 53, No. 8, August 1965, pp. 864- 
876. 

[55] Kouyoumjian, R. G., "The Geometrical Theory of Diffraction 
and Its Applications," Numerical and Asymptotic Techniques 
in Electromagnetics , edTted by R. Mittra, Spring-Verlag, 

New York, 19/6. 

[56] Voltmer, D. , "Diffraction by Doubly Curved Convex Surfaces," 
Ph.D. Dissertation, The Ohio State University, 1970. 

[57] Hong, S., J. Math. Phys. 8, 1223 (1967). 

[58] Levy, B. R. and Keller, J. B., "Diffraction by a Smooth 
Object," Institute of Mathematical Sciences, New York 
University, published under Contract No. AF 1 9(604) -1717; 
December 1957. Also in Comm. Pure Appl. Math., 12, 

1959, pp. 159-209. 

[59] Logan, N. A., "General Research in Diffraction Theory," 

Vol. I, LMSD- 288087; and Vol. II, LMSD-288088, Missiles 
and Space Division, Lockheed Aircraft Corp., December 1959. 

[60] Borovikov, A. A., and Kinber, B. Y., Proc. IEEE, Vol. 62, 

No. 11, November 1974, pp. 1416-1437. 


202 



[61] Fock, V. A., Diffraction, Refraction, and Reflection Waves : 
Thirteen Papers, Air Force Cambridge Research Center Report 
AFCRC-TN-57-102 , June 1975 (AD 117276). 

[62] Pathak, P. H., and Burnside, W. D. , "A Uniform Asymptotic 
Result for the Scattering of Plane Waves by a Circular 
Cylinder," Report 3973-1, March 1976, The Ohio State 
University ElectroScience Laboratory, Department of 
Electrical Engineering; prepared under Contract N62269-74- 
C-0788 for Naval Air Development Center. 

[63] Pathak, P. H., and Burnside, W. D., "A Uniform Asymptotic 
Result for the Scattering by Smooth, Convex Surfaces," 

in preparation. 

[64] Pathak, P. H., "On-Aircraft Antennas," Report 3973-2, 

January 1976, The Ohio State University ElectroScience 
Laboratory, Department of Electrical Engineering; prepared 
under Contract N62269-74-C-0788 for Naval Air Development 
Center. 

[65] Greer, E., personal communications. 

[66] Lipschutz, M. M., Theory and Problems of Differential Geometry , 
Schaum's Outline Series, McGraw-Hill, New York, 1969, 

pg. 196. 


203 



